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ABSTRACT 


This  thesis  will  examine  the  equations  of  motion  for  a 
spinning  ballistic  projectile.  The  goal  of  such  an 
examination  is  to  determine  the  possible  mechanisms  by  which 
a  directed  energy  weapon  may  induce  sufficient  instability  as 
to  significantly  alter  the  projectile's  flight  path.  A 
ballistic  projectile  is  generally  launched  with  a  "fire  and 
forget"  philosophy.  The  desired  impact  point  is  determined 
before  firing.  It  may  be  possible  to  alter  the  projectile  in 
such  a  way  that  it  fails  to  follow  the  desired  trajectory 
thereby  missing  the  intended  target.  Several  variables  appear 
to  be  worthy  of  investigation  to  assess  their  contribution  to 
a  required  instability  or  range  reduction.  Skin  friction  drag 
may  be  increased  from  surface  roughness  generated  by  a  pulsed 
energy  source.  The  results  that  this  thesis  will  examine 
include:  impulse  generated  by  the  laser  interaction, 
additional  Magnus  effects  and  aerodynamic  drag.  Moment 
induced  instability  may  also  result  from  these  in  the  form  of 
a  Magnus  moment  or  drag  torque.  Increasing  the  drag  force 
appears  to  be  the  most  promising  theoretical  solution  to 
defeating  an  incoming  spinning  ballistic  projectile. 
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TABLE  OF  VALUES 


cs 

specific  heat  of  steel 

0.15  cal/gr-C 

d 

155mm  projectile  diameter 

.155  m 

g 

acceleration  due  to  gravity 

9.8  m/sec2 

spin  decelerating  coefficient 

.08 

kd 

drag  force  coefficient 

.004 

kh 

damping  moment  coefficient 

1.2 

kl 

lift  force  coefficient 

1.0 

% 

over  turning  moment  coeff. 

2.5 

kn 

normal  force  coefficient 

1.3 

kt 

Magnus  moment  coefficient 

0.3 

L 

155mm  projectile  length 

70  cm 

lm 

latent  heat  of  melting 

65  cal/gr 

m 

155mm  projectile  mass 

45  kg 

tm 

melting  temperature  of  steel 

1450  C 

u 

wind  speed  over  the  projectile 

600  m/s 

us 

speed  of  sound  in  air 

350  m/s 

V 

speed  of  the  projectile 

600  m/s 

CD 

155mm  projectile  spin  frequency 

220  hz 

Pair 

standard  density  of  air 

3 

1.2  gr /cm 

P steel 

density  of  steel 

3 

8  gr  /cm 

*  Many  of  these  terms  are  approximations  based  on  diagrams  and 
conversations  with  aerodynamics  researchers  at  the  Army  Research 
Laboratory  in  Aberdeen,  Maryland.  The  exactness  of  these  values  is 
not  significant  to  the  magnitude  estimations  in  this  thesis. 
Additional  special  case  values  will  be  introduced  in  the  text. 
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I. 


INTRODUCTION 


This  thesis  is  directed  toward  solving  a  stated  goal  of  the 
Defense  Planning  Guide.  Our  current  armed  forces  are  vulnerable  to 
gunfire  attack  due  to  an  antiquated  counterbattery  system  of 
equipment  and  technique.  Specifically,  in  today's  battlefield, 
surgical  precision  is  required  and  may  involve  protecting  civilian 
lives  in  a  confined  suburban  environment.  There  is  a  need  to 
develop  a  new  system  using  state  of  the  art  technology  to  destroy 
incoming  projectiles  minimizing  loss  of  life.  Laser  technology  has 
given  us  the  ability  to  track,  classify  and  engage  incoming 
ballistic  rounds  to  meet  this  goal.  This  thesis  will  explore  the 
use  of  laser  energy  to  accomplish  the  defeat  of  the  projectile. 

This  mission  responds  to  the  Department  of  Defense  "Defense 
Planning  Guidance  (DPG)  for  FY  1996-2001"  dated  23  May  1994.  This 
involves  countering  armor,  mortar  and  mobile  artillery  threats  to 
U.S.  forces  conducting  an  amphibious  invasion,  U.S.  forces 
defending  a  coastal  region  from  invasion  and  to  coastal  cities  and 
safe  zones  under  U.S.  protection  during  peace  keeping  operations. 
There  are  several  specific  guidelines  in  the  DPG: 

-  "U.S.  and  coalition  forces  would  have  several  key  objectives  in  a  peace 
enforcement  or  smaller  -  scale  operation,  each  of  which  would  require  certain 
capabilities:...  Establishing  and  defending  zones  in  which  civilians  are 
protected  from  external  attacks." 

-  "The  first  priority  in  defending  against  a  large-scale  attack  will  most 
often  be  to  neutralize  the  enemy's  offensive  capabilities." 

"U.S.  military  capabilities  will  be  configured  to  achieve  required 
objectives  while  minimizing  American  casualties." 
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A.  DEFEATING  THE  LAUNCHER  BEFORE  PROJECTILE  IS  A  THREAT 


There  exists  a  need  for  a  Counter -Battery  Weapons  System 
(CBWS)  to  neutralize  the  enemy's  ability  to  launch  large  caliber 
ballistic  projectiles  at  U.S.  forces  and  civilians  under  U.S. 
protection.  The  CBWS  needs  to  be  able  to  determine  the  launchers 
position  from  one  in- coming  round  and  surgically  destroy  that  round 
in  as  rapid  a  manner  as  possible.  A  secondary  purpose  that  this 
thesis  will  not  address  is  eliminating  the  launcher  based  on  cuing 
information  from  the  primary  system. 

The  Combat  Systems  Science  and  Technology  Curriculum  at  the 
Naval  Postgraduate  School  has  a  course  sequence  designed  to  solve 
these  kinds  of  problems  and  propose  system  solutions.  The  class  of 
June  1995,  examined  the  problem  solution  of  defeating  projectile  by 
preventing  the  launch.  The  class  considered  this  problem  to  be 
"shooting  the  archer"  as  opposed  to  "shooting  the  arrow."  A 
complicated  and  expensive  proposal  was  developed  using  acoustic  and 
electromagnetic  sensors  to  detect  an  initial  firing  and  then  a 
guided  projectile  to  destroy  the  launcher  and  nearby  personnel 

B.  DEFEATING  THE  INCOMING  PROJECTILE 

This  weapon  system  should  have  an  integrated,  multi -unit  and 
autonomous  capability  to  designate  the  target  and  automatically 
engage  the  threat  using  a  laser  tracking  beam  and  a  higher  power 
laser  destruction  beam.  Countering  the  threat  from  armor,  mortars 
and  mobile  artillery,  the  system  must  be  able  to  respond  against 
projectiles  ingressing  at  up  to  three  km/ s  from  a  range  of  up  to 
120  km.  The  response  time  for  incoming  projectile  detection 
through  battle  force  threat  warning  and  counterbattery  action  must 
therefore  be  less  than  45  seconds.  The  system  should  be  capable  of 
distinguishing  friendly  fire  from  enemy  fire  through  command 
functions  and  will  also  be  able  to  project  the  most  probable  impact 
point  of  the  in- coming  round  to  determine  if  impact  could  cause 
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major  damage  to  friendly  forces  and  to  provide  warning  of  imminent 
danger  to  ground  forces.  Projectiles  that  pose  no  threat  will  not 
be  engaged.  The  initial  goal  of  the  CBWS  will  be  to  achieve  a  hard 
kill  of  the  projectile. 

Modern  weapon  R&D  is  producing  long  range  over-the-horizon 
(OTH) ,  high  caliber  ballistic  weapons  systems  which  are  small 
enough  to  be  highly  mobile.  These  weapons  systems  are  being 
deployed  in  areas  which  are  becoming  an  important  part  of  the  U.S. 
military  strategy  in  which  Naval  forces  are  being  used  to  support 
ground  forces.  The  shift  from  emphasis  on  open  ocean  warfighting 
to  the  littoral  environment  as  outlined  in  "Forward. .. From  the  Sea" 
will  increasingly  place  forces  at  sea  within  range  of  these  weapons 
systems.  Working  with  the  Marine  Corps  and  their  Operational 
Maneuver  From  the  Sea  (OMFTS)  warfighting  doctrine  will  require 
support  in  various  roles.  These  include:  suppression  of  enemy 
gunfire,  neutralization/denial  of  enemy  operating  areas,  Naval 
Gunfire  Support  (NGFS)  ,  interdiction/neutralization  of  reinforcing 
elements,  and  evacuation  protection.  These  capabilities  will 
require  execution  in  day  or  night,  all  weather,  and  OTH  scenarios 
to  counter  the  threat.  Specifically,  the  threat  must  be  engagable 
in  an  area  lOOnm  wide  and  7  5nm  deep  (from  an  offshore  range  of 
25nm)  .  This  must  be  accomplished  with  an  accuracy  of  several 
meters  to  minimize  collateral  damage  and  fratricide.  The  purpose 
of  the  laser  weapon  system  is  to  provide  the  capability  to  achieve 
these  objectives. 

Within  these  requirements,  our  current  inventory  of  naval 
gunfire  support  (NGFS)  weapons  is  inadequate  to  counter  this 
threat.  The  technological  gains  made  in  other  Navy  systems  has 
occurred  while  naval  guns  have  remained  virtually  unchanged.  The 
maximum  range  of  our  most  capable  system,  the  5"  54 -caliber  M  45, 
falls  far  short  of  the  desired  ranges  and  accuracies. 
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C.  ALTERNATIVE  SYSTEM  SOLUTIONS  TO  DEFEAT  THE  LAUNCHER 


There  are  few  alternatives  to  developing  this  system.  One 
non-material  alternative  may  be  to  adapt  existing  counterbattery 
systems  (Firefinder  radar)  for  use  by  naval  forces.  Because  the 
precision  required  to  surgically  neutralize  enemy  weapon  systems 
does  not  presently  exist  in  the  Navy's  arsenal,  improved  tactics 
cannot  lead  to  a  solution  of  the  stated  problem. 

Currently,  both  the  Army  and  the  Marine  Corps  have  ballistic 
tracking  systems  in  service  to  locate  hostile  guns.  These  again, 
fall  significantly  short  of  the  required  accuracy  and  do  not 
address  the  problem  of  locating  and  destroying  the  incoming 
projectile.  Nevertheless,  these  systems  providing  a  search 
capability  used  in  conjunction  with  a  more  accurate  fire 
control/targeting  system  may  provide  good  initial  or  secondary 
queuing  prior  to  the  terminal  phase  of  the  desired  weapons  system. 

Another  alternative  involves  improvements  to  the  Mk  45 
projectiles.  Several  options  are  being  explored  including 
extending  the  ranges  of  the  projectiles  with  advanced  conventional 
solid  propellants  such  as  liquid  propellant  or  electrothermal 
chemical  propellant.  Further  improvements  to  the  projectile  could 
provide  a  "smart"  seeker  which  would  aid  in  the  terminal  portion  of 
the  ballistic  trajectory.  Use  of  the  existing  TOMAHAWK  cruise 
missile  as  a  weapons  delivery  system  would  require  modifying  the 
missile  to  be  able  to  target  a  point  defined  by  latitude  and 
longitude  vice  using  predetermined  scene  maps  and  an  image  matching 
method.  But,  again,  these  systems  would  do  nothing  to  protect  the 
forces  and  civilians  from  the  projectile  itself. 

D.  THE  SOLUTION  APPROACH  OF  THIS  THESIS 

The  problem  can  then  be  summarized  as  "how  it  is  possible  to 
defeat  an  incoming  ballistic  round?"  The  flowchart  below  describes 
the  major  solution  areas  and  subareas.  This  thesis  is  directed 
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toward  solving  the  problem  of  deflection/destruction  in  flight. 
Specifically,  the  question  "what  can  be  done  to  change  the 
ballistic  projectile's  dynamic  trajectory?"  This  is  accomplished 
by  examining  the  forces  that  control  the  dynamics  and  what  forces 
may  be  applied  to  change  the  dynamics  (Figure  1) . 


Single  Pulse  Multiple  Pulses  Drag  yaw  Torque  Magnus  Effects 


Figure  1  Defining  the  problem  is  the  first  step  in  determining 
possible  solutions  to  defeating  an  incoming  ballistic  projectile. 

The  solution  of  providing  armor  to  defeat  the  round  is 
considered  elementary  and  could  be  considered  as  limiting  the  scope 
of  the  other  solutions  to  larger  caliber  rounds  such  as  the 
155mm.  The  ability  to  defeat  an  incoming  round  in  flight,  and 
ultimately  the  success  of  the  proposed  system,  is  explored  through 
several  theoretical  calculations  involving  laser  interaction  with 
the  surface  of  a  projectile.  Two  main  areas  of  research  appeared 
promising  in  the  problem  solution:  ablative  impulse  and  induced 
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surface  roughness.  The  impulse  problem  was  handily  dismissed  due 
to  the  large  energy  requirements.  However,  the  induced  surface 
roughness  may  have  a  significant  effect  on  a  ballistic  projectile's 
stability  and  range.  This  instability  desired  to  deflect  the  round 
is  most  easily  obtained  by  altering  the  Magnus  moment  of  the 
projectile.  This  instability,  however,  proved  to  be  of 
insufficient  magnitude  to  defend  the  projectile's  target.  The 
additional  drag  created  by  the  increased  surface  roughness  appears 
to  be  a  promising  solution. 

This  thesis  estimates  the  magnitude  of  effects  that  could  de¬ 
stabilize  or  reduce  the  range  of  a  projectile  in  flight.  Most 
previous  work  on  exterior  ballistics  has  focused  on  the  means  for 
achieving  stable  flight  and  maximum  aerodynamic  efficiency.  The 
magnitude  of  the  effects  we  seek  may  demand  great  pointing  accuracy 
and  other  laser  CBWS  problems  may  arise.  These  are  not  addressed 
here . 
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The  155mm  projectile  is  used  as  a  model  to  conduct  these 
calculations  because  it  meets  the  requirements  of  a  medium  to  long 
range  artillery  threat  and  is  a  common  weapon  in  many  international 
forces.  The  155mm  projectile  has  many  variations  of  casing  design 
and  explosive  warheads  that  are  used.  The  most  common  projectile 
is  the  155mm  HE  M107  which  is  used  as  a  general  purpose  round  by 
the  United  States  and  most  other  countries  that  fire  the  155mm 
ballistic  projectile.  The  HE  M107  round  weighs  95  pounds  and  is 
27.5  inches  long.  The  casing  is  composed  of  a  medium  carbon  4340 
steel  (a  carbon,  nickel,  monel  alloy)  that  varies  in  thickness  from 
.25  to  .75  inches  at  the  nose. 


Figure  3  An  exaggerated  description  of  the  increasing  angle  of 
the  projectile  with  respect  to  the  trajectory  along  the 
trajectory. 
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As  the  projectile  travels  through  the  air,  it  spins  at  a 
frequency  of  220-250  hz .  The  spin  is  imparted  by  rifling  within 
the  bore  of  the  launcher .  A  variation  in  the  spin  rate  is  due  to 
the  differing  charges  that  may  be  used  to  fire  the  projectile.  A 
typical  exit  velocity  of  600  meters/second  will  be  used  throughout 
the  calculations.  In  addition  to  spin,  the  ballistic  projectile 
also  undergoes  a  coning  motion  about  the  axis  of  flight  with  an 
origin  near  the  center  of  mass (18.1  inches  from  the  nose  tip) .  The 
angle  of  this  coning  motion  is  typically  3-5  degrees  as  the 
projectile  is  reaching  its  apex  and  then  increasing  somewhat  as  the 
arc  of  trajectory  curves  downward  toward  the  target.  The  coning 
rate  is  between  5  and  20  hz . 

Ultimately,  the  system  can  only  be  successful  if  the 
theoretical  results  can  be  further  refined  and  tested 
experimentally  to  verify  expected  system  performance.  The  purpose 
of  this  writing  is  to  find  possible  theoretical  means  of  using 
laser  energy  in  a  CBWS  to  achieve  the  national  goals  stated  in  the 
DPG .  Specific  details  of  system  design  will  not  be  addressed. 

The  thesis  will  present  arguments  for  five  possible  mechanisms 
of  achieving  a  "kill"  of  the  projectile.  All  of  these  "kills"  will 
be  achieved  by  deflecting  vice  destroying  the  projectile.  The 
first  theoretical  method  discussed  to  defeat  the  projectile  will  be 
momentum  transfer  by  directed  energy  impulse.  A  more  complex 
argument  discusses  the  stability  equation  for  a  spinning  projectile 
and  addresses  the  most  likely  component,  the  Magnus  effect,  in 
relation  to  achieving  instability.  Finally,  increased  aerodynamic 
drag  is  examined  as  the  most  likely  means  to  achieve  a  successful 
engagement  against  an  incoming  spinning  ballistic  projectile. 
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II.  IMPULSE  DELIVERED  BY  LASER  INTERACTION 


A.  LASER  EQUATIONS  AND  IMPULSE  THEORY 


Before  examining  any  calculations  of  ablative  impulse,  it  is 
necessary  to  introduce  the  equations  which  will  be  used.  The 
equations  that  will  be  used  were  developed  in  class  notes  used  at 
the  Naval  Postgraduate  School  for  various  courses  examining 
directed  energy  capabilities  (Schriempf,  1974)  .  The  terms  that 
will  be  used  are  defined  below: 


F Q  =  (l  -  SI)  I  The  laser  power  density  absorbed  where  31  is  the 
Reflectance,  I  is  the  laser  power  density  arriving 
at  the  surface  and  F Q  is  the  power  transmitted 
through  the  surface  per  unit  area. 


tp  The  pulse  length  (duration) . 

D  The  thermal  diffusion  length  is  the  distance 

hr 

required  for  the  temperature  to  drop  l/e  of  the 
surface  value. 


P 


T 


0 


m 


cs,  Cl 
Lm'  Lv 


The  material  density. 

The  initial  temperature. 

The  melting  temperature. 

The  vaporization  temperature. 

The  specific  heat  of  the  solid,  liquid. 
The  latent  heat  of  melting,  vaporization. 
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The  boiling  temperature. 


K  „ 

K  =  —  C 

P 


The  thermal  diffusivity  is  equal  to  the  product  of 
the  thermal  conductivity  and  the  specific  heat 
divided  by  the  density. 


The  criterion  for  vaporization  of  a  diffusion  thickness  is 
given  by : 


Fo  ~n  ~  P  +Lm+cl(Tb~Tm>  +  V  ■ 


(2.1) 


When  numerical  values  are  inserted,  the  Ly  term  dominates  the 
equation  simplifying  to  the  form: 

T, 


(2.2) 


Since,  Dd » 2 /TkO  ,  this  then  becomes: 


Fo-2  “¥=  LvP 


(2.3) 


Solving  this  equation  for  the  pulse  length  gives 

2 

Fv  2 

Tp>4K—  P  . 

Fo 


(2.4) 


Now  to  solve  the  specific  impulse  on  a  target  will  require 
several  assumptions  before  this  can  be  carried  any  further. 

1)  The  beam  radius,  R,  is  much  larger  than  the  diffusion 
length . 

2)  The  target  is  very  thick  compared  to  the  diffusion  length. 

3)  The  melting  will  be  ignored. 
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The  time  for  the  surface  to  reach  vaporization  temperature  is: 


xb  = 


(2.5) 


The  surface  erosion  speed  is  taken  from  equation  2.1 
substituting  A z  for  the  diffusion  length.  This  represents  the 
thickness  of  the  material  evaporated  during  a  time  At.  The 
erosion  speed  is  then  expressed  as: 

'-y  1/[Cs(Tm-T0)-Cp(Tv-Tm)*I.v)  .  <2'6) 


Since  Lm<Lv  ;  Cj«Cs=C,  then  U s *  —  1/  [ Ly+CTV ]  . 


(2.7) 


The  vapor  blow-off  speed  is  derived  from  the  momentum 
conservation,  pvUv  =  ps  Us,  resulting  in: 


pvUv  =  F0/  (LV+CTV)  . 


(2.8) 


The  associated  pressure  with  the  blow-off  is  then: 

2  R  /jj 

P  =  pvUv  =  pv— {gas  constant/molecular  weight)  Tv  so, 

and  combining  this  with  equation  2.8  results  in: 

P  =  F0^]f^/  {LV+CTV)  *  F0JCjT^/[fi{Lv+CTv)]  (2.9) 

since  C«3—  for  most  metals. 

M 

Now,  the  specific  impulse  delivered  during  an  energy  pulse  is: 
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p  xb ^  • 


(2.10) 


Substituting  equations  2.5  and  2.9  into  2.10  and  using  typical 
values  results  in  a  simplified  form: 

XD  ,  2 (2.11) 

Im  =  8£0(l-6-|),  [(taps)  dyn-s/cm  ] 


where  EQ  is  the  pulse  energy  per  unit  area  (fluence) 
=  F0xp,  Tp  units  are  |is . 

This  equation  will  generally  result  in  agreement  with 
experimental  data  within  a  factor  of  two.  This  is  sufficient  for 
the  estimates  of  this  thesis.  If  the  results  of  further 
calculations  indicate  a  possible  effect  within  a  factor  of  two  from 
current  technology,  then  better  estimates  would  be  in  order. 

One  additional  aspect  of  laser  theory  should  be  addressed 
here.  If  the  EQ  is  very  large,  the  delivered  impulse  will  fall  off 
with  increasing  values  of  EQ .  This  is  due  to  the  absorption  of 
laser  energy  within  the  vapor  near  the  surface  of  the  target.  At 
some  value  of  power  density,  the  plume  of  vaporized  material  blow- 
off  from  the  target  will  become  hot  enough  that  the  surrounding  air 
or  vapor  are  sufficiently  ionized  to  absorb  large  amounts  of  the 
laser  radiation.  This  phenomenon  is  called  laser  absorption. 

Another  problem  is  called  beam  decoupling.  When  the  laser 
beam  electron  density  increases  to  a  value  where  the  plasma 
frequency  is  larger  than  the  laser  frequency,  then  the  light 
interaction  with  the  free  electrons  dominates.  The  laser  is  cut¬ 
off  from  the  surface. 

With  these  in  mind,  it  is  important  to  find  the  critical 
energy  density  at  which  these  processes  occur.  The  upper  limit 
density  is  about  1023  electrons  per  cubic  cm  electron  density  where 
the  material  reaches  a  temperature  of  full  ionization.  But,  this 
is  not  accurate  enough  since  the  ionization  degree  is  very 
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sensitive  to  the  temperature.  Figure  4  shows  the  predictions  for 
the  impulse  as  a  function  of  the  pulse  energy  (Nielsen,  1972)  . 
These  values  will  be  used  later  in  the  calculations  for  ablative 
impulse . 


Figure  4  Specific  impulse  delivered  to  solid  targets  by  1.06 
micron  laser  radiation  (Nielsen,  1972) . 
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B.  ABLATIVE  IMPULSE  CALCULATIONS 

l.  Single  Pulse  Impulse 

The  impulse  calculations  are  examined  first  since  this  would 
be  the  easiest  and  most  direct  method  of  deflecting  an  incoming 
round.  A  Laser  induced  ablation  from  a  single  pulse  would  be  the 
preferred  method  of  destruction  if  it  were  physically  possible. 
This  would  require  a  less  sophisticated  system  that  could  "fire  and 
forget."  The  problem  is  complicated  when  a  system  is  required  to 
accurately  track  and  reengage  a  target  multiple  times.  This 
requirement  also  detracts  from  the  ability  to  rapidly  engage  the 
next  subsequent  target  in  a  dense  battlefield.  With  this  in  mind, 
and  using  the  equations  for  critical  energy  of  a  laser  system,  it 
is  possible  to  determine  the  feasibility  of  the  single  pulse 

system. 


14 


To  further  simplify  the  initial  calculation,  a  perpendicular 
impulse  will  be  assumed.  If  the  required  energy  is  at  physically 
attainable  levels,  then  further  calculations  for  reduced  angle 
impulse  will  be  in  order . 

The  impulse  transferred  by  ablation  is  found  by  equation  2.12 
to  be  proportional  to  EQ  for  Tp^lO|is  and  EQ>20J/ cm  when  the  term 
in  the  bracket  becomes  small  compared  to  one: 

Im.8B0(l-(6l|))  . 

E0  (2.12) 


The  impulse  is  calculated  from  the  change  in  momentum  required 
where  the  momentum  of  the  moving  projectile  is, 


p=mv= 1 . 8x10*  gem/ s . 


(2.13) 


For  a  single  pulse  system,  a  shipboard  close-in  defense  would 
require  a  deflection  on  the  order  of  50  meters  for  every  1000 
meters  traveled  by  the  projectile.  The  result  would  be  a  .05  rad 
change  in  direction  or  a  change  in  momentum  of 


(A p)  regd=* 05  ^  ■  8xl°9)  =l08gcin/s. 


(2.14) 


Then,  the  impulse  transferred  in  a  10  square  cm  area  to  the 
projectile  is  estimated  by: 


( Im^  req d' 


(Ap  )  reqd_  i 


=10  dyns/cm‘ 


(2.15) 


Finally,  referring  to  equation  2.12,  the  energy  fluence  required 
will  be  10 °J/cm  with  a  pulse  duration  of  10  6  seconds  Presently, 
this  is  not  within  the  realm  of  possibilities  for  current  laser 
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systems.  Therefore,  further  calculations  of  reduced  angle  impulse 
are  not  required.  However,  there  is  another  way  to  examine  the 
problem.  It  may  not  be  a  close-in  defense  system  that  is  engaging 
the  target  projectile.  Without  placing  any  limitations  on  the 
amount  of  deflection  required  for  a  given  length  of  travel,  a 
calculation  may  be  conducted  for  a  change  of  yaw  angle.  A 
significant  change  of  yaw  angle  of  about  10-20  degrees  would  change 
the  drag  significantly  and  alter  the  range  of  the  projectile. 
Using  a  spinning  top  approach  and  a  small  angle  change,  a  first 
approximation  again  rules  out  the  possibility  of  sufficient  energy 
for  a  single  pulse. 

2 

From  Figure  4,  values  of  Energy  30  J/cm  and  I/E  of  5  dyn  s/J 

2 

show  that  an  impulse  of  150  dyn  s/cm  could  be  attained. 

To  determine  the  required  impulse,  one  finds  the  required 
torque  from  the  required  change  in  angular  momentum  in  the 
following  steps: 


the  torque,  T=Pa, 


(2.16) 


1  2 

the  moment  of  inertia,  I  =  —MR  ,  (2.17) 

and  the  angular  momentum,  L  =  Ito  =  1 0>h?  (2.18) 


Figure  6a  The  torque  applied.  Figure  6b  A  change  in  angle. 


Therefore,  the  torque  can  be  expressed  as: 


dL 

dt 


=  10) 


dn 

dt 


=  1 0) 


da 

dt 


F-a . 


(2.20) 


The  impulse  per  unit  area  required  is: 


^Preg 


F  At 
A 


JO) 

IFI' 


JO) 

aA 


•Aa . 


(2.21) 


Using  the  values  from  Figure  4,  Eczit  =  — with  the  result  of  30 
MW.  This  is  an  energy  fluence  far  beyond  the  laser  absorption 
threshold  and  is  an  implausible  solution  with  current  technology. 

2.  Multiple  Pulse  Impulse 

This  leaves  the  less  desirable  option  of  impulse  delivered  by 
multiple  pulses  to  be  examined.  Using  the  results  of  equation 
2.21,  the  total  energy  required  must  be  divided  into  optimal  energy 
pulses.  If  Eopt  =  30—^,  then  105  pulses  are  required.  Using  10  ps 
pulses  at  optimal  energy,  it  is  possible  to  determine  how  much 
engagement  time  is  possible  as  the  projectile  cones.  The 
engagement  time  is  estimated  by  the  following  using  a  coning  rate 
Of  5  Hz: 


U 

G)c  =  2— Hz  also  v=RcG)c  =  atana  G)c. 


(2.22) 
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With  the  established  coning  rate  and  angular  velocity,  an 
engagement  can  only  be  achieved  within  a  period  of  .06  seconds. 
Applying  these  pulses  over  the  engagement  time  for  a  1  ns  pulse 
results  in  a  power  requirement  of  6xl010  Watts  which  is  still  well 
outside  of  an  expected  capability. 

The  impulse  calculations  demonstrate  that  it  is  currently  not 
possible  to  achieve  a  "kill"  on  a  ballistic  projectile  with  the 
ablation  energy  transfer  alone.  The  next  step  is  to  determine  how 
the  stability  may  be  affected  by  changing  the  aerodynamic  forces 
acting  on  the  projectile. 
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m.  STABILITY 


There  are  several  topics  that  should  be  explored  before 
beginning  the  calculations  to  prove  or  disprove  that  a  viable  CBWS 
can  cause  instability.  Most  importantly,  a  sound  understanding  of 
spinning  projectile  stability  theory  is  established  as  a  framework. 
Then,  using  the  stability  as  a  guideline,  one  term  (the  Magnus 
moment)  in  the  final  stability  equation  appears  to  be  most 
promising  for  further  analysis  in  causing  instability.  The  Magnus 
force  and  ultimately  Magnus  moment  must  be  defined  and  discussed  to 
determine  whether  these  components  are  worthy  of  further 
calculations  to  support  an  instability  argument. 

The  conditions  for  stability  of  a  spinning  ballistic 
projectile  are  given  by  McShane  et  al .  The  following  discussion  is 
a  summary  of  the  results  of  this  analysis  with  the  details  given  in 
Appendix  A. 

A.  THE  EQUATIONS  OF  MOTION 

The  focus  of  this  analysis  is  the  derivation  of  the  equation 
of  motion  of  the  yawing  motion  due  to  the  various  aerodynamic 
forces  acting  on  the  projectile.  The  velocity  vector,  3,  of  the 
projectile  with  respect  to  the  ground  can  be  decomposed  into  a 
component  along  the  projectile  axis  of  symmetry,  u-j_ ,  and  a  complex 
component,  i;  =  u2+iu3,  where  u2  and  u3  are  the  velocity  components 
perpendicular  to  the  projectile  axis  in  a  coordinate  system  fixed 

to  the  projectile.  The  tangential  velocity  of  the  center  of  mass 

F 

is  U.  The  ratio,  is  the  "vector  of  yaw"  which  has  the 

magnitude  of  sin5  .  A  new  variable  is  defined  as: 
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p 

i  fvdp 

X.Ae  °  ,  '3-D 


where  P  is  the  arclength  along  the  trajectory  in  calibers  and  v  is 
the  spin  rate  per  caliber  of  travel.  The  yawing  equation  of  motion 
derived  in  Appendix  A  is  a  second  order  inhomogeneous  differential 
equation  for  the  quantity,  X. 

The  solutions  to  the  homogeneous  equation  are  given  as 
exponential  functions  in  P  as: 


X  = 


M 


If  l  )  ±  (J <2Jr-JA>  ) /o+ivA(l±a)  /BJ  dp 


(3.2) 


where  the  integrand  is  a  function  of  the  relevant  aerodynamic  force 
coefficients.  The  yawing  motion  is  then  given  by  the  superposition 
of  the  two  solutions  of  the  homogeneous  equation  plus  a  particular 
solution  of  the  inhomogeneous  equation: 

X  =  +  L2^2  +  * p •  (3.3) 

The  coefficients  L±,  L2  would  be  determined  by  the  initial  yaw  and 
rate  of  yaw  on  the  trajectory. 

The  projectile  is  stable  along  the  trajectory  if  the  solutions 
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A-j_  and  A2  decrease  in  magnitude  so  that  in  the  limit  of  long 
trajectory,  only  the  particular  solution  remains.  This  means  that 
A  is  independent  of  the  initial  conditions.  The  particular 
solution  need  not  be  small  however .  For  large  spin,  the  rate  of 
precession  and  the  direction  of  its  axis  changes  only  slowly.  As 
the  trajectory  curves,  the  yaw  for  a  stable  projectile  can  still  be 
significant.  In  fact,  the  projectile  becomes  more  stable  along  the 
trajectory  as  the  velocity  of  the  projectile  slows  at  a  greater 
rate  than  the  spin. 

In  reality,  the  yaw  will  be  small  if  the  projectile  is  stable 
and  the  solution  of  the  particular  integral  is  small;  this  is 
almost  always  the  case  unless  high-angle  fire  is  used,  then  the 
particular  solution  may  not  be  small. 

The  analysis  of  stability  of  the  yaw  motion  concentrates  then 
on  the  two  solutions  A1  and  A2  which  can  be  written  in  the  form: 

r  i 

f  [  —  (a1+a4— -)±r]  dp 


where  a1,  a4,  r  are  functions  of  the  aerodynamic  coefficients 
(defined  in  Appendix  A)  and  the  spin  rate  which  is  not  time 
independent  along  the  trajectory. 

One  can  see  that  the  yaw  will  decrease  if  and  only  if  the  real 
part  of  the  first  bracket  under  the  integral  is  negative  and  in 
absolute  value  greater  than  the  real  part  of  the  second. 

In  the  following,  we  will  summarize  this  analysis  with  the 
thought  in  mind  to  determine  which  of  the  aerodynamic  coefficients 
must  be  changed  by  laser  action,  for  instance,  and  by  what 
magnitude  in  order  to  destabilize  the  projectile. 


21 


B.  THE  STABILITY  EQUATIONS 


The  purpose  of  this  section  can  be  achieved  now  by  defining 
the  terms  of  stability  which  are  derived  more  thoroughly  in  the 
appendix.  The  stability  factor  S  and  the  stability  terms  are 
usually  expressed  as: 

1  k 


Hid  ^ 

S1  =  JN+k-2JH~JD~JA~Y~  ' 


ind  ^ 

S ^2  =  rp — — — ,  and 


(3.5) 


S 


3 


2k  2Jh+  (2Jt-2Ja) 


md2 

A 


The  conditions  for 
inequalities : 


stability  are  expressed  as  the 


1 

S 


< 


S2S3 


and  >  0  . 


(3.6) 


It  is  worth  noting  that  the  spin  term,  v,  appears  only  on  one  side 
of  the  inequality.  This  is  useful  for  determining  the  effect  of 
varying  spin  rate  on  a  projectile's  stability. 

These  inequalities  yield  the  desired  framework  for  comments  on 
projectile  stability.  Namely,  there  are  conditions  which  do  not 
satisfy  the  inequalities  that  are  cause  for  instability.  When 
either  of  the  terms  of  equation  3.6  become  negative,  the  conditions 
for  stability  are  violated.  These  are  the  conditions  for 
exploitation  in  the  calculations.  It  becomes  clear  in  analysis 
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(see  the  discussion  at  the  end  of  Appendix  A)  that  the  only  term 
which  can  be  reasonable  exploited  to  cause  instability  without 
catastrophic  damage  to  the  entire  projectile  is  the  Magnus  moment 
term,  JT.  This  is  true  because  this  term  has  a  polarizing  effect 
on  the  first  condition  of  3.6.  Since  this  term  contains  opposite 
signs  in  each  term  of  the  numerator,  then  it  will  have  an  effect  in 
opposite  directions  of  the  over -all  stability.  Further,  because 
the  Magnus  coefficient  can  either  be  positive  or  negative  depending 
on  its  relation  to  the  center  of  mass  location,  it  can  change  signs 
within  each  term  (but,  change  them  both  in  opposite  directions 
simultaneously) . 

The  drag  term  also  appears  capable  of  causing  a  change  in  the 
stability  term,  S2 ,  since  it  is  of  the  same  order  of  magnitude  as 
the  Magnus  term.  However,  the  JN  term  contains  the  drag 

coefficient  as  well.  So  a  change  in  the  drag  will  result  in  a 
corresponding  change  in  JN'  nullifying  any  effect  it  would  have  on 
the  stability.  In  examining  equations  3.5,  it  can  be  seen  that 
increasing  the  drag  will  cancel  out  a  similar  increase  in 
Applying  any  change  in  the  drag  term  to  either  the  denominator  or 
the  numerator  of  equation  3.6  will  result  in  similar  cancelling  of 
the  effect.  This  leaves  the  Magnus  term  as  the  only  term  that  can 
be  altered  in  such  a  way  to  violate  the  stability  criteria 
expressed  in  equation  3.6. 

C.  THE  MAGNUS  FORCE  AND  MAGNUS  MOMENT 

The  Magnus  moment  is  the  most  likely  component  of  stability  to 
be  exploited.  By  increasing  the  Magnus  effect  sufficiently,  one 
can  cause  instability.  Before  making  the  argument  for  instability, 
it  is  appropriate  at  this  point  to  discuss  the  Magnus  force  and 
Magnus  moment . 
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First,  the  coordinate  system  of  the  projectile  should  be 
established.  The  vector  components  of  this  system  will  be  used  to 
define  the  Magnus  terms.  A  ground-based,  right-handed  coordinate 
system  will  be  used  with  a  center  of  mass  of  the  rigid  projectile 
serving  as  the  origin. 


The  projectile  is  considered  a  solid  body  of  revolution. 
Along  the  axis  of  rotational  symmetry,  a  unit  vector  X  is  chosen 
positive  in  the  direction  toward  the  nose  from  the  center  of  mass. 
Therefore  the  total  angular  momentum  of  the  body  can  be  expressed 
as  the  sum  of : 

1)  the  angular  momentum  about  X 
and  the 

2)  total  angular  momentum  about  an  axis  perpendicular  to  X. 
Both  vectors  are  in  the  ground- fixed  coordinate  system. 

The  sum  of  all  the  forces  and  moments  will  determine  the 
position  of  the  projectile.  For  this  discussion,  however,  only  the 
notation  that  is  necessary  to  understand  the  Magnus  contributions 
will  be  introduced  including  several  additional  variables  that 
provide  an  understanding  of  this  effect. 

For  the  angular  momentum  about  X,  the  magnitude  of  the  product 
AN,  where  A  is  the  moment  of  inertia  of  the  body  about  X,  and  N  is 
the  axial  spin  (angular  velocity  with  a  positive  direction  of  a 
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right-handed  screw  about  x  expressed  in  radians  per  second)  result 
in  the  total  angular  momentum  represented  by  the  vector  AN 5c . 

An  aerodynamic  force  or  moment  (torque)  is  produced  by  the 
interaction  of  the  rigid  body  and  the  fluid  of  the  atmosphere 
around  it.  There  are  several  in  addition  to  the  Magnus  components 
that  are  considered  necessary  to  form  a  complete  model  of 
symmetrically  spinning  projectile  behavior.  These  are  outlined  in 
Appendix  A.  The  drag  force  will  be  introduced  in  more  detail  in  a 
later  chapter  as  a  model  for  examining  the  effect  of  a  roughened 
projectile  skin  on  an  aerodynamic  coefficient. 

The  Magnus  effects  on  spinning  projectiles  were  first 
documented  by  G.  Magnus  in  the  1800's.  He  observed  that  a  spinning 
musket  shot  experienced  drift  in  the  direction  of  spin  as  if  there 
were  an  additional  sideforce  acting  on  the  projectile.  It  has  been 
consistently  recognized  that  this  effect  is  an  important 
contribution  in  determining  stability  and  trajectory  of  spinning 
projectiles,  missiles  and  bombs.  Although  there  has  been  a  great 
deal  of  interest  in  these  effects,  there  are  few  studies  that 
examine  these  effects  in  great  detail.  Because  of  the  complexity 
of  this  effect  we  will  at  this  stage  use  the  basic  principle  of  the 
Magnus  force  and  Magnus  moment  as  outlined  below  make  a  simple 
estimate  of  the  effect  of  varying  surface  roughness.  To  date, 
there  have  been  no  models  developed  to  describe  the  complexities  of 
an  asymmetric  flow  caused  by  an  asymmetric  body,  although  one  paper 
addresses  the  non-linear  effects  of  the  Magnus  force  for  large  yaw 
(Platou,  1956)  . 

It  is  widely  accepted  that  the  Magnus  effect  is  produced 
entirely  by  viscous  effects.  An  asymmetric  boundary  layer  (around 
a  symmetric  rigid  body)  is  produced  by  both  the  spin  of  the 
projectile  and  its  angle  of  attack.  Most  of  the  research  to  date 
has  examined  bodies  with  low  angles  of  attack  and  very  smooth 
surfaces.  A  simple  illustration  shows  the  asymmetric  boundary 
layer  and  the  resulting  forces  for  the  low  angle  of  attack. 
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In  Figure  10,  we  consider  a  cross  section  through  the 
projectile  with  its  plane  spanned  by  the  trajectory  and  the 
direction  given  by  £  *  X,  where  €  is  the  vector  along  the 
trajectory.  The  bottom  view  represents  the  projectile  cross 
section  in  this  plane  with  a  boundary  layer  and  the  cross  velocity 
component  of  the  wind  as  it  adds  or  subtracts  to  the  boundary  layer 
velocity  due  to  the  projectile  spin  in  this  plane. 

On  the  upper  portion  of  the  cross  section,  the  fluid  flow  over 
the  projectile  is  parallel  to  the  rotational  motion,  and  on  the 
lower  side  they  are  in  opposite  directions.  On  the  upper  side,  the 
mean  velocity  in  the  boundary  layer  is  larger  than  on  the  lower 
side.  Therefore,  the  pressure  is  larger  on  the  lower  side  and  the 
effect  is  a  net  pressure  force  acting  upward  in  the  Figure.  The 
projectile  is  acted  upon  by  a  net  force  to  the  left  when  viewed 
along  the  trajectory. 

For  a  large  angle  of  attack  (large  yaw  in  excess  of  20 
degrees) ,  the  boundary  layer  begins  to  separate  and  vortices  flow 
into  the  outer  inviscid  flow.  The  result  is  an  equal  and  opposite 
circulation  causing  an  additional  sideforce.  Currently  three 
dimensional  models  of  boundary  layer  behavior  are  too  difficult  for 
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Figure  11  Vortices  shed  due  to  the  Magnus  effect  and  may  lead 
to  additional  force  component  in  the  Magnus  terms  (Power,  1974) . 

theoretical  predictions  beyond  the  crudest  approximations. 

Since  the  effects  on  the  155mm  projectile  are  calculable 
within  a  low  angle  of  attack,  the  further  analysis  of  the  three 
dimensional  boundary  layer  are  not  conducted  in  this  work.  Such 
circulation  and  vortex  shedding  may  contribute  to  rapid  growth  of 
instability  as  precession  occurs  and  should  be  further  investigated 
in  subsequent  work. 

For  the  symmetric  body  before  a  laser  has  induced  surface 
roughness,  the  Magnus  force  can  be  fairly  accurately  described. 
The  Magnus  force  is  the  result  of  the  fluid  interaction  of  the  air 
stream  flowing  over  the  boundary  layer  of  a  yawed  spinning  body. 
The  magnitude  of  the  force  is  pd3JCpNvsin6 ,  where  Kp  is  the 
dimensionless  coefficient  of  the  Magnus  force.  The  force  is 
perpendicular  to  the  yaw  plane.  The  Magnus  force  is  represented  by 
the  vector  for  positive  values  of  N.  The  result  is  illustrated 
in  Figure  12. 


Some  research  has  been  conducted  to  determine  the  value  of  the 
coefficient  and  relates  it  to  the  length  to  diameter  (L/D)  ratio  of 
the  body  for  low  angle  of  attack  and  very  smooth  surfaces  in 
supersonic  flight  (Power,  1974).  Since  we  are  trying  to  determine 
the  effect  of  varying  the  roughness,  this  solution  is  not  helpful 
in  the  end,  but  may  be  used  to  determine  a  ballpark  value  of  the 
unperturbed  coefficient.  The  results  of  this  research  are 
reproduced  in  the  figure  on  the  following  page. 
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Figure  13  Angle  of  attack  vs  Magnus  Force  coefficient 
(Power,  1974) . 

Although  it  is  the  Magnus  force  that  results  from  the 
roughness  that  we  hope  will  ultimately  cause  instability,  it  is  the 
moment  arm  placement  of  this  force  on  the  body  that  causes  the 
torque  about  the  center  of  mass.  This  is  the  contribution  of  the 
Magnus  moment.  If  the  Magnus  force  does  not  pass  directly  through 
the  center  of  mass  of  the  projectile,  then  a  resulting  Magnus 
moment  is  produced.  The  magnitude  of  the  moment  is  given  as  is 
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the  dimensionless  coefficient  of  the  moment.  This  coefficient  is 
approximated  in  McShane,  et  al .  with  a  dependence  upon  the  center 
of  mass  location  in  the  projectile  (the  location  of  the  CM  is  on 
the  horizontal  axis  of  Figure  14)  (McShane,  1953) . 


(McShane,  1954) . 


The  Figure  depicts  firing  data  obtained  for  a  4.5  inch  rocket 
and  although  not  precise  for  our  purposes,  it  provides  a  general 
idea  of  how  the  magnitude  of  the  coefficients  varies.  The  value  of 
the  coefficient  is  read  from  the  graph  at  the  point  where  the 
moments  are  applied  in  relation  to  the  center  of  mass  position) . 
For  the  155mm  projectile,  the  center  of  mass  would  be  near  the 
number  two  on  the  axis  of  Figure  14.  This  gives  a  general  idea  of 


representative  values  for  the  Magnus  moment  if  a  Magnus  force  were 
applied  to  the  forward  third  of  the  projectile  (e.g.  a  Magnus  force 
applied  at  position  4.5  would  yield  a  Magnus  moment  coefficient 
near  zero) . 

The  value  of  the  moment  coefficient  may  be  positive  or 
negative  as  the  mean  position  of  the  force  moves  away  from  or 
toward  the  center  of  mass.  The  hope  is  to  achieve  a  roughness  that 
will  produce  a  force  with  a  different  length  moment  arm  from  the 
center  of  mass  and  thus  produce  a  change  in  the  Magnus  moment 
coefficient.  The  actual  center  of  mass  will  not  have  changed,  but 
the  position  of  the  Magnus  force  relative  to  the  center  of  mass 
will  have  the  effect  of  changing  the  center  of  mass  for  predicting 
the  Magnus  moment  coefficient  with  the  above  Figure.  The  Magnus 
moment  lies  in  the  yaw  plane  and  is  perpendicular  to  x.  The  vector 
[x* (x*v)  ]  has  magnitude  vsin5  .  The  moment  is  depicted  and 
expressed  in  Figure  15. 

It  is  hoped  that  inducing  a  Magnus  moment  by  increasing 
surface  roughness  would  become  a  technologically  feasible 
possibility.  If  it  becomes  possible  to  use  surgical  precision  (to 
within  centimeters)  to  produce  a  roughness  at  a  desired  location 
from  the  center  of  mass,  then  the  resulting  torque  would  be  changed 
by  the  amount  r  X  F,  where  r  is  the  distance  from  the  center  of 
mass  and  F  is  the  magnitude  of  the  Magnus  force.  Ultimately,  then 
the  problem  of  defeating  an  incoming  projectile  could  be  achieved 
with  a  calculated  amount  of  energy  directed  toward  a  specific  point 
on  the  projectile  with  predictable  results.  This  would  be  the  most 
efficient  use  of  the  proposed  system. 
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MAGNUS  MOMENT 


1 

MAGNUS  MOMENT  =  MT  =  pd  KtN[xX(xXv)] 


Figure  15  The  Magnus  moment  (Lieske,  1966) . 

D.  ESTIMATE  OF  REQUIRED  CHANGE  OF  MAGNUS  COEFFICIENT 

The  desired  change  in  the  Magnus  moment  coefficient  can  be 
estimated  by  varying  the  values  of  the  coefficient  and  plotting  the 
results  as  a  function  of  S2,  S3.  This  is  done  in  Figures  16  and 
17  below.  The  first  Figure  shows  how  the  terms  have  a  polarizing 
effect  since  they  contain  opposite  signs  in  S2,  S3.  This  is 
desirable  since  any  change  in  the  value  of  the  coefficient  may 
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cause  an  instability  to  result  from  either  term.  The  second  plot 
is  the  product  of  S2>  S3  demonstrating  the  narrow  region  where 
stability  is  satisfied  (the  region  where  the  product  is  positive) . 
For  typical  values  of  the  coefficient,  a  magnitude  change  by  a 
factor  of  five  may  be  sufficient  to  cause  instability. 


MAGNUS  MOMENT  COEFFICIENT  vs.STABILITY  TERMS  S2  and  S3 


-0.5  0  0.5 


Magnus  Moment  Coefficient 


Figure  16  Stability  factors  vs  Magnus  moment  coefficient  value 


S2  x  S3 


THE  S2  x  S3  PRODUCT 
vs . 

MAGNUS  MOMENT  COEFFICIENT 


Figure  17  The  stability  factor  product  produces  a  region  of 
stability  for  certain  Magnus  moment  coefficients. 

E.  ESTIMATING  YAW  CHANGE  VIA  CHANGES  IN  THE  MAGNUS  MOMENT 

The  next  set  of  calculations  conducted  to  provide  a  theory  for 
instability  are  based  on  a  point  mass  trajectory  model  developed  by 
the  Ballistic  Research  Laboratory.  The  equations  derived  for  this 
model  provide  a  simplified  means  to  calculate  the  yaw  of  repose. 
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ae  ,  a  dimensionless  expression  to  quantify  the  effect  of  changing 
yaw  angle  along  the  trajectory.  The  goal  of  the  BRL  was  to  develop 
a  model  that  would  incorporate  the  effects  of  yaw  without  the 
computing  time  of  a  complex  rigid  body  simulation.  The  estimate 
was  an  intermediate  step  in  achieving  this  goal.  For  this 
discussion,  2e  is  the  desired  term  for  examining  a  magnitude  effect 
of  varying  the  Magnus  moment  (Lieske,  1966) . 

If  the  results  of  this  calculation  demonstrate  a  dependence  of 
the  yaw  of  repose  on  the  Magnus  moment,  then  the  stability  argument 
previously  suggested  is  validated.  The  BRL  found  that  the  cJe  term 
was  successful  in  matching  empirical  data  when  applied  to  equations 
of  motion  that  included  yaw  drag  and  lift.  They  did  not  apply  this 
term  independently  to  varying  conditions  of  the  surface  roughness 
of  the  projectile. 

In  the  development  of  the  yaw  of  repose  estimate,  the  BRL  used 
some  notation  that  is  not  consistent  with  previous  chapters .  To 
avoid  confusion  and  define  some  new  terms,  a  list  of  these  symbols 
is  provided  below  (although  Lieske  et  al .  ,  accounted  for  the 
Coriolis  force,  it  has  been  discounted  by  many  other  studies  as 
negligible  and  is  ignored  in  this  magnitude  analysis  as  well) : 

B  Transverse  moment  of  inertia 

gQ  Acceleration  due  to  gravity  at  ground  level 

T  Unit  vector  in  the  direction  of  v 

Kn  Drag  force  coefficient 

o 

Kn  Yaw  drag  force  coefficient 

ua 

The  same  coordinate  system  previously  introduced,  a  ground- 
fixed  right  handed  Cartesian  system,  will  be  used.  The  body  is 
still  considered  a  solid  of  revolution  and  the  unit  vector  x  is 
along  the  axis  of  symmetry  pointing  forward  through  the  nose.  The 
magnitude  of  the  angular  velocity,  N,  is  parallel  to  51  and  is 
positive  for  a  right-handed  screw  motion  along  x. 

The  following  are  the  simultaneous  differential  equations  of 
motion  for  a  spin  stabilized  projectile.  First,  for  the  center  of 
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mass : 


Pd2 

m 


(K  +&2Kd  ) 

^  o  CC 


V*  ( xx  v) 


pd\,  -r  ,  Pd3^ 

- — K  vx  + - 

m  s  m 


(x*v)  +g. 


(3.7) 


Recall  that  the  total  angular  momentum  H  can  be  expressed  by: 

H  =  ANx+B  (x*£)  .  (3.8) 

Therefore,  the  sum  of  the  applied  moments  gives  the  vector 
rate  of  change  of  the  angular  momentum: 


5  =ANx+ANX+B  (x*X)  =p  d3KMv(v*x)  -pd^Wx^+pd^^x*  (x*v]  -p  diKANvx. 


(3.9) 


Defining  the  yaw  of  repose  from  the  rigid  body  system  of 
differential  equations  requires  these  assumptions: 

1)  The  projectile  can  be  represented  sufficiently  as  a 
body  of  revolution. 

2)  The  projectile  is  dynamically  stable. 

3)  The  initial  yaw  is  assumed  to  be  small  and  has  a 
negligible  effect  on  the  trajectory. 

In  equations  3.7  and  3.9  the  magnitude  JX*  17 1  is  present  in  the 
lift  and  Magnus  terms.  If  we  let  I=— ,  then  the  magnitude 
Ixxj|=sin5.  Lieske,  et  al.,  then  define  the  yaw  of  repose: 

3esfx  (xx i)  =x-Cos5 J.  (3.10) 

This  vector  represents  the  vector  yaw  directed  from  T  toward  J?  and 
is  in  the  plane  determined  by  I,  x.  The  effect  on  the  trajectory 
by  ae  will  be  small  with  the  assumed  conditions.  Therefore,  (Ee 
is  assumed  negligible  yielding  the  following  approximations: 
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5  =  cos  8i,  5  =  cos8j.  (3.11) 

Then,  dividing  the  total  angular  momentum  into  parallel  and 
perpendicular  components  results  in: 

AN=-pdiKfNv,  ANX+B(x*X)  =  pd3JCMv(  v*x)  -pd4KHv(x«5)  +pdiKTN[x*  (x*v)  ]  . 


Then  substituting  x  and  its  derivatives  in  3 . 7  and  3.12  with 
equations  3.10  and  3.11  will  provide  a  solution: 


il= 


-pd2  {KD+b2KD)  v2I  pd2KLy2ae  pd3Ksvcosbl  p d2KpNv 


m 


m 


m 


m 


5/1)  +£T . 


(3.13) 


ANcosfi  J +Bco s 6  (a9x?)  +Bcos28  (1*1)  =pd3JCMv2  +pd4JCrWv[cos8  (Ss+cos 6l)  -f]  -pd4Ktfvrcos6  [  (a^+cosCJ)  *2]  . 


Cross  multiplying  both  sides  of  equations  3.13  with  I  and  utilizing 
the  terms,  I,  I  the  yaw  of  repose  becomes: 


[ -AKjNcosft  (?*7)  +ind2K'rNcos6  [?*  {u-§+pdiK3cos6^r/m)  ]  ~KL[v-  {v'l)  T]  [2Bcos35  {vT)  -pd*KHcos26v3]  ] 
‘  [pd3KIiC^v4  +  pd5J<'yi!C3jW3cos5v3+iCLBcos6  [  (v*f)3-  (v-v)  ]  ] 


(3.14) 


With  the  v  term  dominating  the  denominator,  the  last  term  can  be 

2 

discarded,  leaving  a  ratio  of  .  The  middle  term  in  the 

v4 

denominator  after  Kt  can  be  considered  approximately  equal  to  g. 
With  some  known  behavior  of  artillery  shells  this  further  reduces 
to  a  final  form: 

2  > 

-AKjJsKvxti)  +md  KyN[v*  {ti-g)  ]  (3.15) 

2  = -  . 

0  O  yl  r-  /"X 

Pd  KLKMv  +  P d  KpKjN  v 


The  next  step  is  to  use  this  equation  to  estimate  the  effects 
of  varying  the  K  terms.  v*^  is  estimated  using  the  ideal 
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2  .  2  * 

trajectory,  this  is  done  by  approximating  ae=sm  0 
can  be  written: 


Then ,  3.15 


ae  =  (ab-aa)  (v* u)  -a b(v*§)  , 


(3.16: 


such  that,  Q.e*(tta-2ab)  (vQ*§) 


(3.11] 


Now  the  useful  form  desired  for  sin  5  can  be  written: 


2  *  2  2,2  2  . 
sin  8  =  vQg  ( a a+4ub-4 a aab) 


13.18) 


For  the  mathematical  analysis  it  is  necessary  to  define  the 


terms  <xa,  ocb  . 


a*  =  — 3 


akln 


;3.i9) 


p  d  3  KlKmv 4  +  p  ^KpK^v2 


iaKrjH 

^  jb  .  o  o 

pdJCL%v  +pd  iCpK-yN  v 


It  is  worth  noting  that  again,  both  terms  contain  coefficients 
for  the  Magnus  moment.  Using  numerical  values,  and  changing  the 
Magnus  moment  coefficient  from  0.3  to  1.5,  the  solutions  for  the 
equations  3.19  are: 


aa  «  9x10 


-9  sec' 


a*.  *  4x10 


-8  sec' 


Thus  the  result  for  initial  values  of  yaw  will  be: 


<2a  ~  ixio  5  *  sin25,  8  ~  /10’xio  2. 

Similarly,  if  the  Magnus  moment  is  increased  by  a  factor  of  five, 
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then  the  results  are: 


aa  ~  9»10  9 ,  ab  ~  7xio  7 ,  Sg  *  lxio  3,  b  ~  /TOxio  1. 

This  is  an  order  of  magnitude  change!  The  theoretical  evidence  to 
support  instability  with  a  five-fold  magnitude  increase  in  the 
Magnus  moment  is  established.  These  results  although  not  linked 
directly  to  the  stability  criteria  support  the  argument  that 
instability  can  be  achieved  with  a  change  in  the  Magnus  moment. 
Based  on  the  stability  equations,  we  expected  instability  to  arise 
by  changing  the  Magnus  moment.  Here,  we  have  changed  the  Magnus 
moment  and  verified  that  the  yaw  angle  will  change  considerably. 
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F. 


WILL  THE  MAGNUS  EFFECT  CHANGE  SUFFICIENTLY  WITH 
INCREASING  THE  SURFACE  ROUGHNESS 


A  drawing  of  the  projectile  as  viewed  from  the  nose  tip  coming 
out  of  the  page  will  be  helpful  for  estimating  the  magnitude  of  the 
Magnus  effect. 


The  cross  flow,  v± ,  is  approximately  50  m/s  for  a  small  yaw  of 
about  five  degrees.  This  is  nearly  equal  to  the  velocity  of  the 
boundary  layer  close  to  the  projectile  surface  (laminar  flow)  which 
is  calculated  to  be  60  m/s.  Applying  the  Bernoulli  equation  to  the 
unroughened  laminar  flow  projectile  surface  leads  to  the  following 
analysis : 


=  constant  =  PQ  + 


(3.20) 


where  v-^v^  giving 
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4] 


Adding  a  surface  roughness,  as  in  Figure  19,  to  the  entire 
body  would  cause  the  resultant  velocity  vector  on  the  lower  side  to 
decrease  as  fluid  momentum  from  outside  the  boundary  layer  is  drawn 
into  the  boundary  layer  by  turbulence.  The  velocity  along  the  top 
of  the  projectile  will  not  change  appreciably  since  the  component 
of  cross  wind  velocity  and  the  velocity  of  the  boundary  layer 
(caused  by  the  angular  velocity  of  the  spinning  projectile)  are 
nearly  equal.  The  result  is  a  small  change  in  the  Magnus  force 
that  will  be  estimated  now. 

For  the  projectile  with  a  symmetric  band  of  roughened  surface, 
the  view  from  the  nose  shows  how  the  turbulence  draws  the 
surrounding  flow  into  the  boundary  layer  in  the  roughened  region. 
The  turbulence  transports  fluid  momentum  across  the  boundary  layer 
and  tends  to  bring  the  flow  velocity  closer  to  the  surface.  Notice 
the  turbulence  does  not  change  the  profile  of  the  velocity  along 
the  top  surface  of  the  projectile. 

The  turbulent  boundary  layer  is  turbulent  at  a  Reynolds 
number , 

„  . p  v- 


in  a  band  around  the  projectile,  using  the  full  flow  velocity 
including  the  longitudinal  component  because  it  is  the  full 
velocity  at  the  boundary  layer  edge  that  determines  the  Reynolds 
number . 
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unchanged  profile  on  top 


Then  comparing  the  equations  in  3.20  results  in  the  following: 


P2  +  |pv2  =  P0  +  |p  vl 


(3.21) 


approximating  then 


1.21-2.  _  132 

p2  *  Po*  2P<’'x-4Vi>  *  POt2P'4VJ. 


such  that  now, 

A pturb  “  ■gPvi. '  "*  ^p  laminar~^ptuzb  ~  ~~q  Pvx  ’ 

This  is  a  change  from  the  laminar  case  by  a  reduction  of  about  25% 
change  in  pressure  and  hence  a  similar  change  in  force.  But,  this 
is  only  applied  to  the  band  of  roughness  which  for  a  laser  induced 
case,  may  only  be  one  tenth  of  the  length  of  the  projectile.  As  a 
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moment  arm,  this  would  represent  20%  of  the  moment  arm  length. 
Thus,  applying  a  25%  decrease  in  Magnus  force  to  only  20%  of  the 
surface  would  result  in  about  a  5%  change  in  the  Magnus  moment. 
This  is  not  a  significant  enough  effect  to  achieve  the  instability 
we  had  hoped  for  via  the  Magnus  moment. 

If  the  roughness  were  applied  to  the  entire  surface,  the  net 
result  would  only  be  a  change  in  the  Magnus  force  by  a  total 
magnitude  of  about  25%.  This  is  still  far  short  of  the  five-fold 
increase  which  is  likely  to  cause  an  instability.  There  may  be 
additional  non-linear  factors  that  contribute  to  the  Magnus  effect 
such  as  those  described  by  the  vortices  that  produce  an  additional 
sideforce  (Platou,  1956) .  The  analysis  applying  this  additional 
component  may  still  provide  validity  to  the  prospect  of  exploiting 
the  Magnus  moment  as  a  means  to  achieve  instability,  but  is  not 
undertaken  here . 

The  case  of  an  asymmetric  roughness  does  not  present  a  uniform 
force  that  could  be  applied  to  the  Magnus  effect  analysis.  In  fact, 
the  asymmetry  may  only  cause  a  wobble  of  the  same  proportions 
derived  above.  The  asymmetry  may  also  appear  symmetrical  to  the 
Magnus  effect  if  the  spinning  boundary  layer  creates  a  uniform 
turbulence  at  high  spin  rates.  However,  it  may  create  a  moment 
arm  of  drag  which  will  be  discussed  later. 
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IV.  AERODYNAMIC  DRAG  AS  A  MODEL 


The  drag  will  be  examined  as  the  friction  factor  changes  by  a 
factor  of  five;  this  value  was  not  chosen  arbitrarily,  however  the 
analysis  of  the  Magnus  effects  proved  that  a  similar  change  in  drag 
would  not  result  in  an  instability.  As  such,  it  is  hoped  that  a 
similar  magnitude  change  in  the  drag  will  result  in  an  additional 
drag  force  sufficient  to  reduce  the  range  of  the  projectile.  Drag 
will  be  examined  separately  to  determine  if  the  effects  of  drag 
alone  can  prevent  an  incoming  ballistic  projectile  from  reaching 
its  target. 

The  model  for  the  analysis  will  be  a  cylinder  to  approximate 
a  projectile.  The  flow  through  a  cylinder  has  been,  in  particular, 
very  well  studied.  A  pipe  is  a  cylinder  with  internal  flow  that 
approximates  the  flow  behavior  on  the  external  surface  as  well.  It 
has  been  well  established  that  increasing  surface  roughness  in  a 
turbulent  flow  will  have  an  effect  on  friction  resistance.  In  this 
case,  for  the  projectile,  the  effect  is  on  the  aerodynamic  drag 
across  the  whole  body.  For  the  model,  then,  we  want  to  increase 
friction  resistance  by  a  substantial  factor  of  three  to  five. 

The  friction  resistance  is  due  to  the  surface  effect  of 
changing  the  boundary  layer.  By  increasing  surface  roughness,  the 
sublayer  will  break  up  and  become  increasingly  turbulent. 

For  the  155mm  ballistic  projectile  a  roughness  value,  e,  of 
.046mm  (White,  1979)  can  be  divided  into  the  length  (700  mm  is  used 
as  an  estimate  to  account  for  the  geometry)  to  enter  Figure  21  for 

friction  with  smooth  and  rough  walls.  With  a  value  —  of  15000 

7  e 

and  using  a  typical  Reynolds  number,  10  for  the  projectile,  the 
flow  can  be  determined  as  past  the  transition  point,  but  not  fully 
turbulent.  Therefore,  the  condition  of  turbulence  is  met  and  an 
increase  in  surface  roughness  should  result  in  dramatic  changes  in 
the  drag.  In  the  next  section,  the  calculations  will  be  made  that 
show  the  actual  values  of  the  Reynolds  number  and  how  the  roughness 
must  vary  to  provide  the  necessary  results. 
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Figure  21  Reynolds  number  vs  drag  coefficient 
(White,  1979) . 


ft..  AERODYNAMIC  DRAG  ESTIMATES 


The  Reynolds  number  for  our  projectile  must  first  be 


estimated : 


Bk  «2.6*107 


(4.1) 


■L  V 


Then,  entering  Figure  21  with  that  result,  and  moving  acr< 


the  chart  with  the  value  of 


=  15000  to  determine  the  drag 


coefficient,  gives  CD  =  .004.  An  increase  of  more  than  three  can 
be  seen  on  this  Figure  using  CD  £  .013  resulting  in  an  of 
approximately  200.  This  results  in  a  required  roughness  of  3.5  mm, 


e  =  ttvv  =  3  ■ 5mm  (4.2) 

200 

Similarly,  if  the  drag  coefficient  is  used  as  the  measure  to 
determine  a  change  in  roughness  required,  it  is  possible  to 
calculate  the  coefficient  by  (White,  1979)  : 


KD  = 


■  031 
R  1/7 


8700 


(4.3) 


Using  this  result  and  multiplying  it  by  five  will  again  result  in 
a  change  of  roughness  L/e  on  the  order  of  several  hundred  and  a 
resulting  roughness  of  several  millimeters.  This  roughness  can  be 
visualized  as  dimples  covering  the  surface  of  the  projectile  that 
are  several  millimeters  deep  and  spaced  so  that  no  smooth  surface 
is  larger  than  the  diameter  of  those  dimples.  The  number  of 
dimples  required  can  be  calculated  using  the  area  of  the  projectile 
and  a  dimple  size  of  5  square  centimeters  to  be  approximately 
162,000.  This  means  that  the  laser  must  pulse  that  many  times  on 
the  incoming  projectile.  With  a  short  pulse  length  on  the  order  of 
nanoseconds,  this  can  be  done  in  milliseconds.  A  longer  pulse 
length  but  no  longer  than  microseconds  will  still  allow  sufficient 
roughening  in  under  a  second.  With  the  coning  and  spinning  motion 
of  the  projectile,  a  pulse  rate  somewhere  in  the  middle  will  allow 
for  the  most  dispersion  of  pulses  over  the  entire  body. 

Similarly,  a  continuous  pulse  laser  may  melt  the  surface 
sufficiently  in  less  than  1/2  second  to  cause  the  same  effect.  The 
roughness  does  not  need  to  be  caused  by  dimples.  It  may  be  equally 
effective  to  produce  ridges  of  the  same  magnitude  of  solidified 
molten  flow  produced  by  the  laser  energy  for  a  roughening  effect. 

From  the  criteria  for  vaporization  equation,  an  approximation 
of  the  energy  required  can  be  made  (Schriempf,  1974) : 
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F0Xp  ~  P&x(Cs{Tm  Tq)  +lM'  F0Xp  ~  (56  00  -),  (4.4) 

cm 

where  Ax  is  the  hole  size  or  roughness  desired.  Substituting  in 
the  Ax  =  0.35  cm  term  results  in  the  following  requirement 
(details  are  located  in  Appendix  B) : 


The  symmetry  of  the  roughness  does  not  significantly  affect 
the  net  force  acting  on  the  projectile.  Because  the  projectile  is 
spinning  and  the  roughness  will  be  asymmetric  over  the  entire  body, 
the  turbulent  boundary  layer  will  spin  with  the  projectile  as  the 
cross  wind  interacts  with  the  turbulence.  In  fact,  an  asymmetric 
disturbance  of  several  thousand  holes  around  the  circumference  of 
the  spinning  projectile  will  produce  the  largest  increase  in  the 
drag  effect  since  the  boundary  layer  and  hence  surface  friction 
will  increase  around  the  entire  surface. 

B.  DRAG  VS.  RANGE  ESTIMATE 

The  question  remains,  "how  does  this  drag  reduce  the  range  of 
the  projectile?"  To  complete  this  estimate,  the  equations  of 
motion  incorporating  the  drag  will  be  used.  The  vertical  drag  will 
be  assumed  negligible  compared  to  the  axial  drag  for  ease  of 
estimation.  This  is  probably  a  reasonable  assumption  anyway  since 
the  projectile  will  be  engaged  close  to  its  terminal  point  and  the 
altitude  will  be  small  compared  to  the  range.  Additionally,  the 
yaw  as  the  projectile  reaches  its  target  may  be  at  its  maxima 
providing  additional  lift  to  counter  any  vertical  drag.  This 
assumption  also  allows  a  starting  calculation  for  the  time  of 
projectile  flight  which  will  be  needed  in  the  more  detailed 
equation  of  motion. 
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v y  t-  — gt  . 


(4.5) 


Using  a  firing  angle  of  45  degrees  for  the  projectile  to  achieve 
maximum  range,  the  horizontal  and  vertical  components  of  the 
velocity  are  equal  initially  ( Vv  =VY  )  .  The  projectile  is 
calculated  to  have  79.5  seconds  as  a  time  of  flight. 

The  equation  of  motion -used  to  determine  range  information  is: 


x  =  - - p  ( Y) - j-r-xKp  =  -G(v)  x. 

m  u5(Yj  u 


The  functions  that,  in  general,  are  dependent  upon  altitude  will 
also  be  considered  constant  for  the  estimate  made  here.  Applying 
the  following  steps  to  equation  4.6  results  in  a  final  form  useful 
for  determining  the  range  shortfall  with  an  increase  in  drag: 


Finally, 


so  that 


x  =  u,  u  =  -u G(v)  ,  In  U\t)  =  -g(v)  t, 

u(o) 


u(t>  =  e-u(o)G(v)t.x(t)  =j-e-G(v)  tdttXo_ 


-G(v)  t. 


x(  t) 


V  _  V 

xo  -G{v)  t  „  U  T,  xo 
- e  +X„,  where  X  =- - 

G(v)  °  °  G(v) 


o  .  ^  -G(v)  t. 


x(t)  =  —If  (1  -  e 
G(v) 


Applying  the  final  form  of  4.7  to  the  projectile's  flight  in 
two  steps  will  yield  the  desired  results.  First,  by  choosing  a 
point  2000  meters  from  the  CBWS  the  projectile's  time  of  flight  and 
range  up  to  that  point  can  be  calculated.  Then,  using  that  point 
as  the  start  of  the  engagement  with  a  the  new  drag  coefficient,  the 
range  of  the  engagement  leg  of  flight  can  be  compared  to  the  range 
if  the  drag  were  not  changed.  The  result,  for  our  increase  in  drag 
is  a  20%  reduction  in  range  with  the  engagement  starting  at  2000  m. 

From  the  launch  of  the  projectile  until  2000  meters  away  from 
the  target: 
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For  the  remaining  2000  meters  of  travel  new  velocities  in  both 
the  vertical  and  horizontal  components  must  be  applied  in  addition 
to  an  increased  drag  coefficient  (increased  from  .004  to  .02) .  The 
velocity  will  be  calculated  from: 

x(t)  =  vx  e  t  (4.8) 


where,  G=  .002283  for  a  drag  factor  of  .004.  This  will  give  an 
additional  projectile  range  during  the  increased  drag  portion 
(G=  .008121)  of  its  flight  path  of  only  1619  meters.  The 
projectile  will  fall  short  of  the  target  by  nearly  400  meters  or  20 
percent  of  the  engagement  distance.  A  twenty  percent  reduction  in 
range  with  an  easily  obtainable  laser  requirement  is  an  encouraging 
result ! 

C.  YAW  TORQUE  CREATED  BY  A  CHANGE  IN  ROUGHNESS 

In  the  case  where  only  one  pulse  or  a  series  of  pulses  in 
close  proximity  roughen  the  surface  of  a  projectile  in  a  small  area 
asymmetrically,  there  will  be  a  torque  induced  that  will  produce  an 
increase  in  yaw  with  the  same  frequency  as  the  spinning  projectile. 
Depending  on  the  existing  yaw  angle  and  the  cone  rate,  this  torque 
may  cause  tricyclic  motion  of  the  projectile. 

If  the  asymmetry  is  out  of  phase  with  the  cone  rate,  it  will 
cause  a  coning  motion  within  the  larger  coning  motion,  thus  a 
tricyclic  motion  as  in  Figure  22.  This  will  almost  always  be  the 
case.  There  will  be  points  in  the  tricyclic  motion  where  the  yaws 
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are  complimentary  and  result  in  a  larger  yaw  and  points  where  the 
yaw  motions  cancel.  The  effects  of  this  complex  motion  are 
difficult  to  describe  and  would  vary  in  almost  every  case. 

With  this  in  mind,  the  projectile  will  be  assumed  to  have  zero 
cone  for  this  calculation  initially.  This  will  enable  a 
calculation  of  the  contribution  of  the  roughness  independent  of  the 
other  forces. 

Dividing  the  drag  force  into  components,  the  component  along 
the  axis  will  cause  a  yaw  to  develop.  The  torque  created  by  this 
force  is  solved  by  the  following  equations: 

T  =  a*F  =  131  lFlsin<l>  (4.9) 
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=  I  51  I Fl sin (180-P) 

*  151 lFlsin(180-Y) 

=  I  51  I  Fl siny 

=  I  a  I  I  Fl  — 
a 

«  r*F,  where  r=7  cm  «  J? . 


Figure  23  Asymmetric  drag. 


An  equivalent  expression  is  then  used  to  solve  the  equation  in 
terms  of  flight  time  remaining  to  the  target.  At.  The  change  in 
angular  momentum  is  related  to  a  change  of  yaw  angle  that  is  chosen 
large  enough  to  increase  drag  sufficiently  to  miss  the  target.  The 
result  is  expressed  below  as  a  sequence  of  steps  assuming  there  are 
three  seconds  of  flight  time  remaining: 


F-r 


Al  .  Ls0)Att 
At  At 


where  Ls  = 


(4.10) 


This  results  in  a  required  force  of  approximately  300  Nt.  But, 
this  is  not  a  useful  number  unless  we  can  attain  it  with  our 
estimated  roughening  that  increased  drag  by  a  factor  of  five.  For 

this, 

F  =  —  pU2CD-A,  SO  A F  =  pu 2  (A Cd)A  *  35  Nt.  (4.11) 

2  u  * 


This  result  is  an  order  of  magnitude  less  than  the  300  Nt  required. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  stated  goal  of  attaining  a  system  that  will  perform  the 
requirements  set  forth  in  the  Defense  Planning  Guide  is  within  the 
realm  of  physical  possibilities.  The  use  of  directed  energy  may  be 
a  viable  solution  to  the  problem  of  defeating  an  incoming  spinning 
ballistic  projectile. 

The  following  potential  mechanisms  were  investigated  with 
available  theoretical  models  and  quantitative  estimates  were  made 
with  the  following  results: 

(1)  Impulse  transfer  of  momentum  of  10%  of  the  projectile 
momentum  by  laser  ablation:  The  allowable  laser  energy  fluences 
before  laser  absorption  waves  will  decouple  the  laser  from  the 
target  are  orders  of  magnitude  too  small  to  achieve  the  required 
lateral  momentum  change. 

(2)  Impulse  angular  momentum  change  by  laser  impulse  to 
increase  the  yaw  angle  significantly  to  induce  large  drag:  Because 
of  the  limitation  of  laser  pulse  energy  by  the  absorption  wave 
threshold,  it  would  require  multiple  sub-nano  second  multi -gigawatt 
pulses  to  achieve  significant  angular  momentum  changes. 

(3)  The  stability  of  the  projectile  and  its  dependence  on 
various  aerodynamic  coefficients  was  investigated  using  a  well 
known  stability  model.  The  Magnus  force  and  Magnus  moment  were 
identified  as  the  most  likely  candidates  for  affecting 
destabilization  by  a  five  fold  increase  of  the  Magnus  force.  It 
was  shown  that  such  an  increase  would  indeed  lead  to  a  significant 
change  in  the  yaw.  However,  a  simple  analysis  of  the  effect  of 
surface  roughening  by  a  factor  of  3-5  by  laser  melting  or  dimpling 
leads  only  to  a  25%  change  in  the  Magnus  force  because  the 
turbulence  tends  to  wash  out  the  pressure  differential  between  the 
top  and  bottom  of  the  projectile. 
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(4)  The  effect  of  increasing  aerodynamic  drag  by  laser 
roughening  of  the  projectile  surface  was  investigated.  It  was 
found  that  a  3-5  fold  increase  of  surface  roughness  (drag)  would 
require  fluences  of  2-5  Kilojoules  per  square  centimeter  and  that 
the  effect  of  the  increased  drag  could  lead  to  a  20%  decrease  in 
the  projectile  range  over  the  last  2000  meters. 

(5)  Creation  of  a  yaw  increasing  torque  by  a  single 
asymmetric  high  roughness  laser  melting  spot  on  the  projectile  was 
investigated  and  it  was  found  that  with  a  fivefold  increase  of 
roughness  on  a  10  square  centimeter  spot,  the  torque  is  too  small 
to  cause  significant  yaw  increase  and  subsequent  drag  increase  and 
range  reduction  in  the  time  available  for  engagement. 

The  aerodynamic  drag  was  certain  to  increase  regardless  of  the 
Magnus  effects  when  the  surface  of  the  projectile  was  roughened. 
The  result  although  not  dramatic  emphasizes  the  viability  of  using 
a  directed  energy  CBWS  as  a  means  to  achieve  a  "kill"  of  an 
incoming  enemy  projectile.  If  the  target  of  the  projectile  were  a 
ship,  then  the  range  reduction  would  certainly  be  sufficient.  For 
land  targets,  the  decrease  in  range  may  pose  other  problems  in  a 
congested  or  civil  warfighting  environment. 

The  system  requirements  may  be  of  a  scope  that  is  entirely 
different  in  nature.  It  is  necessary  to  develop  a  system  that  is 
capable  of  tracking  a  projectile  for  ranges  of  several  hundred 
meters  during  an  engagement.  This  must  be  done  with  an  accuracy  of 
centimeters.  The  system  must  also  be  capable  of  determining  the 
exact  location  of  the  projectile  relative  to  the  projectile's 
intended  target.  It  would  not  serve  its  purpose  to  engage  a 
ballistic  projectile  that  was  certain  to  overshoot  its  aim  point  if 
doing  so  would  cause  it  to  fall  short  and  on  target.  These  are 
just  some  of  the  system  problems  that  will  require  addressing 
before  a  suitable  CBWS  could  be  built. 

The  energy  requirements  for  the  laser  to  induce  sufficient 
drag  are  well  within  the  bounds  of  current  technology.  In  that 
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light,  it  is  advisable  to  find  an  experimental  means  of  validating 
these  findings.  There  may  be  some  interest  in  determining  exactly 
how  deep  the  laser  effects  will  penetrate  or  what  size  melt  ridges 
will  develop  on  the  ballistic  projectile  surface  when  atmospheric 
disturbances  have  reduced  the  effectiveness  of  the  laser.  A  sea 
environment  in  particular,  may  have  a  deleterious  effect  on  the 
directed  energy  against  a  relatively  low  flying  target  such  as  a 
terminal  ballistic  projectile. 

Other  recommendations  for  research  include  a  more  detailed 
look  into  the  drag  behavior  on  a  projectile.  The  model  of  a  flat 
cylinder  does  not  take  into  account  such  things  as  the  point  where 
the  shock  wave  detaches  from  the  body.  There  may  be  other  unique 
aerodynamic  behaviors  with  various  projectile  body  shapes  that 
would  require  significantly  more  or  less  roughening  to  achieve  the 
desired  results. 

Finally,  one  last  attempt  to  exploit  the  Magnus  effects  is 
worth  the  effort.  The  nonlinear  effects  of  the  vortices  shedding 
may  be  sufficient  to  produce  the  desired  instability.  It  may  be 
more  desirable  to  pursue  the  use  of  a  kinetic  energy  vehicle  to 
cause  a  change  in  the  center  of  mass  location  such  that  the  Magnus 
moment  causes  instability.  This  would  likely  be  a  solution  that 
would  provide  a  much  larger  safety  margin  of  destruction  from  the 
projectile's  intended  target. 
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APPENDIX  A 


Aerodynamic  Forces  Acting  on  the  Projectile  Preliminaries 


CP  s  center  of  pressure 
R  =  resulting  aerodynamic  force 


Figure  24  A  single  resultant 
force  on  the  projectile. 


For  this  introduction  into  the  stability  a  preliminary  model 
will  be  used  to  define  some  of  the  most  basic  motions  and  forces 
acting  on  the  projectile.  The  projectile  in  the  Figure  above  is 
rotationally  symmetric  with  velocity,  u  ,  with  respect  to  the  air 
and  an  angular  velocity,  0)  ,  equal  to  zero.  The  terms  of  this 
Figure  are  defined  below: 

M  =  Overturning  moment  =  component  of  T  along  the  Z-axis 

D  =  Axial  Drag  =  component  of  ft  along  the  X-axis 

L  =  Lift  Force  =  component  of  5  along  the  Y-axis 

M,  D,  and  L  depend  on  the  yaw,  6,  by: 

L=X sin6  X,  m  are  independent  of  6 

M=msin8cos8 

D=pu2d2KD,  L=pu2d2XLsin8  ,  M=pu2d3XMsin8cos8 
X=pu2d2KL,  m=pu2d3KM 

D,  m,  and  x  are  all  functions  of  the  projectile  shape, 
diameter,  velocity,  sound  speed  and  density. 


a 


Figure  25  The  forces  and 
coordinates . 
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Since  it  is  desirable  to  have  the  projectile  travel  with  some 
yaw,  it  is  necessary  to  impart  a  stabilizing  spin.  This  will  lead 
to  a  complete  description  of  the  forces  acting  on  the  projectile  by 
reducing  the  aerodynamic  forces  to  a  finite  number  of  coefficients. 
These  have  been  verified  experimentally  within  a  satisfactory 
degree  of  accuracy  for  describing  the  motion  of  the  projectile  with 
a  small  yaw. 


?((3,u)=  net  aerodynamic  forces  acting  at  the  center  of  mass, 
3  ((S ,  u)  =  net  torque  due  to  aerodynamic  moments. 

,%2  unit  vectors  fixed  to  the  projectile. 


F=Fixi 

G=Gi;£i 

U=Uj_^£ 
CO  =(])±g 


CM. 


Then  solving  for  F £  and  G^  is  done  with  a  Taylor  expansion 
about  the  unyawed  position.  The  coefficients  of  the  expansion  are 
independent  of  u  and  (0  to  the  first  order.  Adding  s,  the  angle 
of  rotational  symmetry  between  zero  and  pi,  the  complex  numbers, 
£=u2+iu3  ;  T)  =o>2  +i^3  where  u2=— (5+0  i  u3=— become 
necessary.  Linear  functions  of  u2 ,  u3 ,  (02 ,  (i>3  can  be  written  as 
linear  functions  of  £,£,T|,T).  The  Taylor  expansion  then  becomes: 


F1  =  al  +bi%  +  +  ^2^ 

^  =  F2+iF3  =  a2  +c 15  +  c2i 1  +  d +  d2  rf 


(A.l) 


The  coefficients  except  a  are  in  general  complex  and  the  special 
form  of  F  is  necessary  because  it  must  be  real  for  all  values  of 
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£,r |  .  Rotating  the  vectors  ,&2  >^3  through  the  angle  s  about  the 
axis  of  symmetry  changes  the  axis  system  to  ,£ 2  '^3  anc^  the 

following  transformations  take  place,  £  -  I;  e  1S  ,  T|  -  r|  e  1S  ,Z  Z  e  1S 
rf  -  r|  e  ,  F1  -  F1 ,  and  &  e  .  Because  of  the  symmetry,  such 
a  rotation  cannot  change  the  coefficients  in  equations  A.l.  The 
resulting  equations  are: 


al  + 


e"is  + 


b2r\ 


els  + 


b±le 


is 


T,  1S 

b2r\e 


(A. 2) 


^  e  1S  =  a2+c1i  eis  +  c2x]  e  +  d^e  is  +  d2rie  is 


Comparing  equations  A.l  and  A. 2,  the  terms  on  the  right  from  before 
and  after  the  transformation  can  be  set  equal  in  F-j_  for  all  T)  . 
Choose  r|=0,  £=  h1  -  2b1  b1=2b1  I^cos s,  coss^l  -  b-^0,  similarly  b2=  0. 
Using  the  same  equations  to  compare  & ,  the  d  terms  can  be  examined 
the  same  way  and  they  result  in  coefficient  values  of  zero  as  well. 
Similar  arguments  can  be  made  for  the  components  of  G.  The  result 
is  that  if  a  projectile  is  symmetric  about  its  axis  of  rotation, 
its  force  and  momentum  system  can  be  written  to  a  first  order 
approximation : 

Fl=ai  ir=F2+iF3  =  C-^+C 2r\ 


(A. 3) 


^=G2+iG3  =  C3£+C4n 


where  a1,e-^,c1  ...c^  are  functions  of  density,  shape,  velocity, 
angular  velocity,  and  wind  speed  and  are  independent  of  the 
orientation  of  the  axes  ,St2  ,£3  with  respect  to  the  body.  One 
can  further  show  that  if  the  projectile  has  a  plane  of  mirror 
symmetry,  then  in  A.  3  the  terms  F± ,  Stc-L ,  9tc4  ,  %c2 ,  Sc3  are  even 
functions  of  and  ^  ,  9tc2 , 9lc3 , $c-|_ ,  Sc4  are  odd  functions. 
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The  resulting  dimensionless  form  of  the  force  and  moment 
coefficients  are: 

C- l  =  -pd2u1KN+  ipd3co1iCF 

4  .3 

C*2  =  p  C?  G)  1  ^ JKF  P  ^ 

=:— pC?  G)^iC(ji"“-ipC/  ( A-  •  4 ) 

4  .  5 

“  p  c?  u1-PCH  +  -i  p  c/ 

*1  =  "Pd^i Kda 

Gi  =  -pd4UlKA. 

The  K' s  are  all  even  functions  of  &)-,  and  of  — ,  d—  .  The  signs 

1  a  u 

are  chosen  so  that  for  a  typical  projectile  all  coefficients  would 
be  positive  (for  most  shells,  KT  is  negative)  . 

The  forces  and  moments  are  then  summarized  as : 

2 

pd  u i-Km£  The  Normal  Force  cross  force  due  to  cross 

velocity. 

3 

ipd  ( £  The  Magnus  Force  Magnus  cross  force  due  to 

cross  velocity. 


pd 

ipd3u1F5ti 
p  d  CO  K  j,£ 
ipd3u1K^, 
pd 


The  Magnus  Moment 
Over  Turning  Moment 


Magnus  cross  force  due  to 
cross  spin. 

Cross  force  due  to  cross 
spin. 

Magnus  cross  torque  due  to 
cross  velocity, 
cross  torque  due  to  cross 
velocity. 

cross  torque  due  to  cross 
spin 


Damping  Moment 
or  Yawing  Moment 


5 

ipd  yrrp 

,2  2 

pd  u1KDA 

4 

pd  U-LtO-LiC^ 


Magnus  cross  torque  due  to 
cross  spin. 

Axial  Drag 

Spin  Decelerating  Moment 


The  force  R  has  been  resolved  into  forces  along  and 

perpendicular  to  the  projectile  instead  of  u.  Taking  into  account 

2 

the  yaw,  5,  the  normal  force  becomes  pd  (ucosfi)  X^(usin8)  ,  and  the 
axial  drag  becomes  pd  (ucos5)  KDA .  Computing  these  with 
components  along  and  perpendicular  to  the  velocity  vector  gives: 


pdzuzKLsin6=pdzuzKNsin5cos  5-pd  u  JTDAsin8cos*8 


pd2u2KD=pd2u2KNsin25cosb  +  pd2u2KDAcos^6  . 


(a. 5; 


A  comparison  then  gives: 


2  2 
K ^=K j^COS  8  -Kjjj^COS  8 


(a. 6; 


2  3 

K 6cosb+KDAcos  8 


which  to  a  first  order  approximation  ( 5<1 )  is: 

KL~KN~KDA  an<^  KD~KDA'  (A.l) 

While  the  moment  term,  KM,  remains  identical  to  the  one  previously 
introduced. 

The  general  equation  of  motion  for  a  symmetrical  projectile 
has  four  assumptions : 
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1)  A  fixed  earth  three  dimensional  coordinate  system  ( X ,  Y,  Z) 

neglects  the  Coriolis  force. 

2)  Values  such  as  density  that  change  with  altitude  are 

assumed  constant. 

3)  There  is  no  wind  velocity. 

4)  The  projectile  has  an  angle,  s,  between  zero  and  pi  of 

rotational  symmetry  and  a  plane  of  mirror  symmetry. 

The  velocity  vector,  angular  velocity  vector  and  projectile 
fixed  coordinate  system  remain  the  same  for  the  notation  of  this 
section.  Equations  A. 3  represent  the  Aerodynamic  force  and  the 
Aerodynamic  torque,  respectively.  A  unit  vector  parallel  to  the 
vertical  axis  in  the  earth  fixed  system  is  defined  to  be 
Yj_  Y=Y'2+^1’3  and  the  acceleration  due  to  gravity  transforms 
~gY1  -  -gT .  The  remaining  terms  to  introduce  are:  the  moment  of 
inertia  about  the  axis  of  the  projectile.  A;  the  moment  of  inertia 
about  any  axis  perpendicular  to  the  projectile  axis  through  the 
center  of  mass,  B;  and  the  mass  of  the  projectile,  m. 

The  general  equation  of  motion  is : 

^(m,u)  =F-mg9.  (A.8) 

The  angular  momentum  equation  in  the  projectile  coordinate 
system  is  given  by: 

2  =AC013?1  +BG)2^2  +  5(0 3-^3  =  (A-B)  0)1^1+B(5 

(A. 9) 

and  —7-  =  G  ,  -j-{{A-B)  +BG>}  =  G. 

dt  dt  11 

The  time  derivative  of  vector  P  is  expressed  in  the  projectile 
system  where  P  transforms,  P-  {P-^;  P  =  P2+iP3) .  Using  index 

— ♦  y\ 

notation,  P=P^^  and  P=P j£j+P jX .  The  vectors  £j  are  fixed  in 

the  projectile. 


P  =  PjJtj+Pj  c5 
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Xj  =  (x>x£j 


(A. 10) 


l*k  -  tjOjk'Pj  -*k  -  pk~pj 


In  components,  the  first  one  becomes: 


P-£ 1  =P1+  [  ( P2+1p3 )  («2'<03)  -(P2-iP3)  (w2  +  i(03)  ^  ~o  =  ^1+ 


(A. 11) 

The  second  and  third  added  together  with  the  third  multiplied  by  i 
becomes : 

P^2  +  iP‘;?3  =  P2+iP3~iPl  (°>2  +  i0)3^  +itol  <P2+iP3^  =  P-iPi^l  +i(01P  (A.i2) 


So,  the  equation  of  motion  in  projectile  based  coordinates  are: 
m[u1-(^ri-fTi)  -|]  =F1-mg1 


i[$-iu1T|+i<D15]  =&'-mgT  =  C1^+C2ri -mgT. 


(a. 13; 


And,  the  equations  for  the  angular  momentum  are: 


(A-B)  co1^1+B(3  =G. 


(a. 14; 


The  above  lema  on  the  time  derivative  gives 


—  [(O-L^]  ={  (0*fi-0*Tl )  —  =  Cb-L 

{  0  +  0  =  -IW-lTI 
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d  d. 

—  [o5]  ={  <i)1+(rrri-ri-ti)  -  =  61 

{  tf|  +i(0T|  =  tf|  . 

With  this  we  have  the  two  equations: 

( A~B )  C«3  +  B(jd-^  = 

(A. 15) 

(A-B)  (-io^)  +  B-1\  =  &  =  C3Z  +  C4T1 


The  differential  equation  for  a  change  in  the  force  of  gravity  in 
the  projectile  frame  is  given  by  assuming  Y  =  constant ,  Y  =  0  in  the 
lab  frame  and  applying  the  lema: 

-£z[Y]  =  {  Y±  +  (Tti-TT|)4  (A.  16) 

at  x  2 

{  T  -  iY 1r|  +  icoT. 

There  is  however  an  altitude  dependence  on  the  speed  of  sound 
and  the  density.  The  altitude  change  of  the  projectile  must  be 
tracked  in  order  to  insert  the  correct  p(Yj,  a  (Y)  values.  If  u  has 
a  component  ft?  in  altitude  direction  then  its  altitude  change  is 
given  by: 

Y  =  u-Y  =  U-l^  +  u2Y2  +  u3Y3  =  u1Y1  +  ~  (5T  +  IT)  .  (A.  17) 

The  normal  equation  of  motion  is  then  derived.  Again,  several 
simplifications  are  introduced  to  get  an  approximate  trajectory 
before  making  corrections: 

1)  The  projectile  axis  is  tangent  to  the  trajectory. 

2)  The  only  forces  acting  on  the  projectile  are  drag  and 

gravity. 

3)  The  variables  affected  by  altitude  are  considered  constant. 

4)  There  is  no  wind. 

The  projectile  will  remain  in  the  vertical  plane  of  the  initial 
velocity  vector  as  indicated  by  Figure  26. 
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Figure  26  The  drag  opposes  the  motion. 


D  =  pd2U2KD,  u2  =  X2  +  Y2,  D  -  pd  U  Fg  = 

The  equations  of  motion  are  then: 

mi  =  -pD  uKdu  -  mgf  (A.  18) 

with  components: 

wit  =  -pd ZuKdX,  mt  =•  -pd2uKDt  -  mg  (A.  19) 

where  KD  =  ,  u2  =  X2*?2,  P=P  (r)  ,  a-a(Y>  .  Most  exterior 

ballistics  involves  the  precise  integration  of  equations  A. 19. 

The  standard  form  of  a  normal  solution  is  given  by. 

X  =  ucos9 ,  t  =  usin0 

and  X  =  iicos0  -  usin0  0  ,  Y  =  iisin0  +  ucos0  0  ,  so  that 

mucos0  -  musin0  0  =  -pd^XnU2cos0  (A. 20) 

2  2. 

and  miisinO  +  mucos0  0  =  ~pd  XpU  sin0  -  mg. 

Multiplying  the  third  and  fourth  equations  above  by  cos  and  sin 
respectively,  then  adding  the  result  gives: 


(A. 21) 


2  2 

mu  =  -pd  KDu  -  mgsin0. 

Similarly,  multiplying  by  sin  and  cos  respectively  then  subtracting 
gives : 


0  =  -  — COS0 
u 


(A. 22) 


To  the  above  equations  A. 21  and  A. 22  the  following  terms  are 
added:  p  =  p(y)  ,  a=a(Y)  ,  KD=KD(~)  . 

<3. 

Next,  the  angular  motion  of  the  projectile  will  be  discussed 
starting  with  the  simplification  of  the  equation  of  motion, 
equations  A. 13  and  the  following: 

Acb  ^  =  G-^ 

Bi 1  +  ( B-A )  iO)1T|  =  % 

Y±  =  (Tri  -  ft])  ±  (A.  23) 

T  =  iY1n  -  ico-i^T 
3>  =  *  |  (IT  *  IT) 

where  SP  =  C ^  +  C2T) ,  =  C35  +  C^v\  ,  and  the  values  C^.  .  .C^,  F1 ,  G-j_ 

in  terms  of  the  aerodynamic  coefficients  are  given. 

•,  Changing  the  notation  somewhat  to  include  a  density  factor, 

p - ,  is  necessary.  The  density  factor  is  applied  to ,  the 

m  .  .  dJ 

aerodynamic  coefficient  to  result  in  J  terms,  J-.  =  p - Kn. 

m  cP 

Further,  we  will  define  the  spin  per  caliber  travel  as :  v  =  <0-, - . 

U1 

The  coefficients  then  become: 

U1 

C±  -  (  +  ivJpjm-^j- 

C^  =  (  V  J ^  1  J^) 

O3  =  ( —v ‘ iji  ~  iJ^jmu^  (A.  24) 

“  (—  J ^ijt)  mu-^d 

G-^  =  ~J m  • 

Then  use  approximation  by  means  of  the  normal  equation  of 
motion.  Let  the  center  of  mass  motion  be  governed  by  the  normal 
equation  of  motion: 
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-pd2u2-^  -  gsin0  =  -u2-^  -  gsin0 


-geos- 


(A. 25] 


X  =  ucos0 ,  Y  =  usin0, 


The  solution  is  assumed  to  be  known  and  satisfies 


U  =  -JD— j-  -  gsin0, 


(A. 26) 


Next,  we  will  conduct  a  change  of  independent  variables 
starting  in  the  yaw  equations.  First  replace: 


r  U1  r  U 

p’f  ~ddt-f{  d)dt- 


(A. 27) 


This  is  then  the  arc  length  along  the  trajectory  in  calibers.  The 
equations  determining  the  yaw  motion  will  become  essentially 
independent  of  the  size  of  projectile.  A  large  projectile  has  the 
same  period  of  yaw  measured  in  calibers  of  travel  as  a  small  model. 


(f»  =  =-  Q'  — 

dP  d  dP  d  u  d 


(A. 28) 


The  trajectory  normal  equations  are  then: 

0  0 

u'  =  -JDu  ~  gdsin-y.,  0'  =  -gdcos — 

U  jjZ 


(A. 29) 


The  equation  for  the  derivative  of  the  angular  velocity  becomes 


Ati)  2_  =  Gq  W  2. 


dG- 1  J  3  TT^  TJ  U  V 

—  =  -dvm—  4-  with  o,  -  -^v,  (O'  =  V'-  +  -U' 
AU  A  U  1  d  d  d 
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u  V  U  V 

Therefore,  v  —  +  — -C7'  =  - dxmJ A —  ,  v  '  =  -  —  [/' 
n  H  ri  2i  rr 


2  ^A 

d  vm-^-  resulting 


v'  =  ^ md 2  +  vjd  +  vgdd  2sin0. 


(A. 30) 


To  change  dependent  variables  the  variables,  5,  r| ,  T  are 
changed  as  follows : 


A.  I.  0 


L  J  vdp 


J. 

ij’vdp 


y  =  gTe 


d  0 


gi  =  g*i 


/vdp 


—  is  the  vector  yaw  which  has  the  magnitude  of  the  sine  of  the 
angle  of  yaw.  The  factor 

P  t 

e[ijvdp]  =e[iJco1dt] 

0  0 

changes  the  coordinate  system  from  ^  /  £3  which  rotates  with 

the  projectile  to  one  that  does  not  rotate  about  the  jS^-axis. 

The  quantity,  A  ,  is  the  complex  yaw  measured  from  the  axis  of 
the  projectile  to  the  trajectory  on  a  coordinate  system  with  one 
axis  along  the  projectile  axis  and  which  is  not  spinning  about  the 
axis.  Similarly,  |i,  y,  g,  have  interpretations.  The  coordinate 
system  thus- described  is  determined  when  we  fix  the  direction  of 

^2  • 

The  aerodynamic  force  perpendicular  to  the  projectile  axis  is, 

l 

*m  {-JN  +  i  v  J mt7—  +  ( v  J J-J ^ )  mUr\ 


and  in  the  non-rotating  system  the  lateral  force  is  this  quantity 
multiplied  by 


69 


p 

ifvdp 

0 


This  equation  then  becomes 


=  (~JN  +  ivJF')I nlj2~k  +  (vJXF  +  lJs')mu2^  ■ 


(A. 31] 


The  component  of  drag  in  the  plane  perpendicular  to  the 

projectile  is  due  to  an  aerodynamic  force  acting  on  the  projectile 

( 1 ) 

with  the  complex  number  representation,  &  ,  together  with  that 

component  of  the  axial  drag  which  is  in  the  '  Plane.  This 

2  iT2 

component  is  \JnU  —  .  The  aerodynamic  force  perpendicular  to  the 

u  d 

projectile  has  the  complex  number  Representation  described 
by:  [  ( jp  -  +  ivjp)k  +  (vJxf  +  T^e  real  axis  is 

perpendicular  to  the  trajectory  pointing  downward  and  the  imaginary 
axis  is  horizontal  pointing  to  the  left.  The  quantities  A  and  (i 
are  the  cross  velocities  divided  by  the  axial  velocity  U  and  the 
cross  angular  velocity  multiplied  by  d/U. 

The  final  form  of  the  equation  is  then  obtained.  For  the  sake 
of  abbreviation,  the  term  p 

fvdp  =  <P  •  Making  the  same  transforma¬ 
tions  for  i; ,  1*1,  T  :  0 

1=1'^.  =  (Ac7e~i<p)  'P  =  [A'U+A  (-Jn-gdU~2sinQ)  u-iyku]  e'1(p-^ 


=  [A  '-ivA-A  ( Jn+gdU  sin0)  ]  — — e 


P 


(A. 32] 


Similarly, 


f|  =  [|j/-iYH-|i( JD+gdC7-2sin0)  ] — -e  2<p 

a 


(A. 33) 


and  gT  =  (y'-ivy)-^e  2<^ . 


(A. 34) 


These  are  now  inserted  into  equations  A. 23  leading  to  the  following 
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equations : 


m[X  ' -X  ( JD+gdU  2sin0)  -i^]  =  {C1'i+C2‘i\  ~mgT)  =C1X  —  +C2|i—  -iffy— 5- 


U  U  ‘  n2 


2  d  d 

B[\i'-\i(JD+gdU  sin0)  ]  -iAv(i=C3A  — +C4n— 


(A. 35) 


y'=ig\i,  g1'=(y\i-y\i)  - 


B  2 

Substituting  in  the  constants  C-^.  .  .C^,  Glt  - -sic  (the  transverse 

radius  of  gyration  in  calibers)  gives:  m<^ 

X  —  { J J 3.v J gdU  sin0)  A  +  g+i)  [A— ydu 


_2  —2  —2  —2  v 

\i'=  {-vJT~iJM)  Xk  +(JD~k  JH+ivk  JXT+gdU  sin0+iA— )|i 


Y'=isrH,  Sri'=(Yki-Y^) 


(A. 36; 


Finally,  the  solution  of  the  equation  of  motion  will  be 
completed  by  using  some  approximations  which  will  simplify 
equations  A. 36: 

1)  Each  J-term  contains  a  density  factor  that  is  very  small  in 

-4 

magnitude,  on  the  order  of  10 

2)  A/B  will  be  approximately  1/10. 

ad  -3 

3)  As  with  the  density  factor,  -^-<10 

4)  Derivatives  of  the  J- terms  ar^  negligible. 

The  equations  A. 36  will  be  abbreviated  as  follows: 

X '  =  a1X+a2\i+b 

I x'  =  a3A+a4ix  (A. 37) 

with  the  constants  equal  to: 

al  =JD~JN+1VJE 
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a2~VjXF+i 

(  — V  J rj-t—  hj^)  k 

a4~  (JD~k~2JH+gdU  2sin 0+i^v] 
i>=  (gAsin0-y)  du  2  . 


Next,  we  will  eliminate  H  from  equations  A. 37  by  substituting 
the  first  into  the  second  equation: 


A  "+A'(-a1-a4- 


l2  a2 
—  ')  +A  (a-]_a4 -a2a3 +a-]_ — 


,  .  Z  ,  jO 

ai  ,+i<^r  +a4'i 


To  deal  with  the  derivative  term  a2  ' ,  recall  equation  A. 30.  The 
a2  '  2 

terms  —  ,  a-^ '  are  of  the  order  of  J  ,  whereas  the  other  terms  in 
the  coefficients  are  of  the  order  of  J.  So,  they  can  be  neglected. 
The  result  is: 


A  '  '+A  '  (-a1-a4)  +A  {a1a4~a2a.2)  +b(a4~j-  ')  =0 


(A. 38) 


The  solution  of  the  homogeneous  part  of  A. 38  is  done  with  a 
change  of  variable  to  eliminate  the  first  derivative  term  (the 
particular  solution  will  be  addressed  later) : 


\f  (a1+a4)dp 


A  -  ge 


(A. 39) 


From  this  equation,  the  two  derivatives  A  '=  [q'+-^  (a4+a4)  q]  e a  and 
k"=[q"+ (a±+a4)  q'+^- (a1+a4)2q+^  (a±+a4)  'q]  ea  are  found. 
Substituting  these  into  A. 38  gives: 


q"-i2q  =  0 


(A. 40 : 


where  r 

=  —  [  (a1+a4)  2 -4  (a1a4-a2a3)  ~2  (a-|_+a4)  ]  =—  [  (a]_~a4)  +4a2a3~2  (a-^+a4)  ] 
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To  find  the  approximate  solution  of  equation  A. 40,  let 

Z=  (lnq)  '  =  —  q';  Z'  =  --^-q'2  +  —' '  =  -Z2+r2.  Then,  Z'+Z2-r2=0.  If  r 
Of  q  o 

were  constant  then  tne  equation  would  have  a  solution  Z=±r.  r  is 

not  really  constant  because  of  the  v  in  the  a1...a4  terms,  but  r 

is  nearly  constant  which  allows  a  trial  solution  of  Z=r+e,e<r.  It 


follows  e  '+2re+e  +r  '=0  . 


r  '  2 

Let  6=-— —  then  neglecting  6.' ,  e 


quantities  satisfies  the  equation  approximately.  The  error  is  of 
2  2 

the  order  (e '+e  )/r  .  The  solution  is  then: 


Z~r~- 


and  hence, 


Using  the  other  sign  in  Z=±r , 


/t-s? 


(A. 41! 


(A. 42) 


g  =  e 


Solving  for  A ,  we  also  expect  two  solutions 


j[a1+a4-—  )/2±r]dp 


(A. 43) 


A  =  e 


Now,  to  express  the  coefficients  in  A. 43  in  terms  of  aerodynamic 

coefficients,  the  coefficients  ^2.+a4  ^rom  A. 37  can  be  summed  and 
.  2 

also  substituted  into  the  value  of  4r  found  in  A. 40.  The  result 
is  aided  by  using  A. 30  and  yields: 


4  r  -  A  -+4k  JM+  [JN  Jp- 
Ba 


■k  2Jtt-{2 JT-J 


.  .  md‘ 


2  iAV 


(A. 44) 


Several  specialty  terms  are  required  to  complete  the 
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computation  of 


The  stability  factor, 


A2v2 


4BjC  Jr 


0N1_5 


as  well  as 


(A. 45] 


„  ,222 
E(4r2)  =  A  . 

s2 

r  ' 

Neglecting  the  imaginary  part,  if  —  is  computed  precisely,  then 

the  denominator  multiplied  by  its  conjugate  complex  with  a 

magnitude  comparison  leads  to  the  same  result  as  if  neglecting  the 

2 

imaginary  part  of  4r  (except  for  r«l  which  is  a  special  case  not 
discussed) . 


(4r2> (vV) du-2sin9*-h. 

r  2  ( 4 r 2 )  2v2o2  v  °  °  A  A 


(A. 46: 


The  yawing  motion  is  then  a  combination  of  the  two  solutions 


fe  ° 


+iv4)±('4iv2t4Jt'2v(^-JD-*'V(2Jr-JA)  ^ >  -tt>  2idp 


(A. 47) 


An  approximation  for  the  solution  of  a  stable  projectile  is 
such  that  the  ratio  of  the  imaginary  and  real  parts  of  r— <—  for 

£  U 

stable  projectiles.  Using  the  binomial  theorem  to  express  the 
square  root, 


'  =a+b  -  2r=i  —  va+  [JN~JD~k  ZJH~  (2JT~JA) — —  ]/a. 

ID  /i. 


(A. 48: 


In  that  case,  near  stability  the  yaw  motion  is  given  by  a 
combination  of  the  two  solutions: 
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0  i/[  <J»-JD-k'2-  <Z-V*V  +  (1±0)  1  dP 


(A. 49; 


This  was  the  form  that  was  used  at  Aberdeen  Proving  Grounds  for 
reduction  of  data.  An  alternate  form  for  a  stable  projectile  can 
be  written  as  follows: 


/■  ,  2 .  /  ,  2  2  2  2  1 2 

r  (4r  )  A  vv  a  v  a  .  .  T  T  md  .  ,  2 

—  = - 2~“- - T-(" - 5 —)  (JD-JA—T^a  ' 

1  2(4r  )  BZ  BZ  A 


(A.50: 


and  then  using  equations  A. 43,  A. 44,  and  A. 48  leads  to:  (A. 51) 

■1  j [2  JD-Js-ic"2  Jj/+gdcr"2sin0-  ( JD-md2  +gdC7"29in0)  / o2±  [ J„- JD-k '2 JH-2  JT-  JA)  ]  /o+  [ivA (l±o)  B]  ]  dp 


The  criteria  for  stability  can  now  be  adequately  discussed.  A 
projectile  is  deemed  to  be  stable  if  small  disturbances  have  no 
permanent  effect,  that  is,  if  the  yawing  motion  does  not,  in  the 
limit,  depend  on  the  initial  yaw  and  the  initial  yawing  motion. 
The  yawing  motion  of  a  projectile  is  given  as  a  linear  combination 
of  the  two  solutions  of  A. 43  plus  a  particular  solution,  so  that 

X  “1,^  X  2.  "*■  1^2  X  ^  X  p  . 

The  projectile  is  stable  if  X^_,  A2  decrease  in  magnitude  so 
that  in  the  limit  only  the  particular  solution  X ^  remains  (it 
becomes  independent  of  the  other  terms) .  The  particular  solution 
need  not  be-small  though,  for  a  large  spin,  the  rate  of  precession 
is  small,  and  the  direction  of  its  axis  will  change  slowly.  Since 
the  trajectory  curves  downward  after  the  apex  of  flight,  the  yaw 
will  become  larger.  The  question  then  becomes,  "what  are  the 
necessary  and  sufficient  conditions  that  the  solutions  given  by 
A. 47  approach  zero  as  the  arc  length,  p,  increases?"  One  can  see 
that  the  yaw  will  decrease  if  and  only  if  the  real  part  of  the 
first  bracket  under  the  integral  is  negative  and  in  absolute  value 
greater  than  the  real  part  of  the  second. 
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Let 


a-(  JD+JN+k  jh  ' 

A2V2  -2 
b=~  —  ^  +4ic  2Jm, 


(A. 52: 


A  ”2  /nd 

c=2  —  v  [JN-JD-k  JH-  (2Jt-Ja)  — —  ]  so  that 


q  2  J 

k:  ^ e  0 


J  [  (-a+iv-^)  ±  {b+ic)  2  ]  dp 


The  condition  of  stability  is: 


a>  IMy/b+icI  . 


2  2 

Reformulating  the  stability  criterion  b+ic=\Jb^ +c  (cos<J>+isin(j>) 


where  cosd)=—  ^  and  sind)=— — — 


From  De  Moivre's  theorem: 


^2+c2 


b2+c2 


>+Icl  =4^b2+c2  -i<1  + - = - )  =*  (b+\/^2+c2)  ■ 

\ 2  ,/i^  N  2 


(A. 53: 


So  that  the  stability  criteria  become: 


a>*  (b+Jb^c2) 


where 


-a=R (a1+a4) 


b=R4r ' 


c=Jm  4r  . 


Additional  steps  are  taken  in  reformulating  the  stability 
criteria  in  terms  of  aerodynamic  coefficients: 
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and  finally. 


>  {b+sjb2+c2)  , 

a>0 

2a2  >  b+]jb2  +c2 , 

a>0 

2  a2-b  >  \jb2  +  c2 , 

a>  0 

{2a2 -b)  2  >  b2  +  c2 

,  a 

4  2.2 

4a  -4a  b  >  c  , 

a>0 

2  2  4 

-4a  b  >  c  -4a  , 

a>o 

2  c2 

~b  >  -a  2  +  -— — , 

4a 

a>0 . 

(A. 54) 


2  2  c2  2 

Noting  that  a  is  of  the  order  of  J  and  - a  is  of  the  order 

4 

,2  _  t2  2  ,.2 

A  2  u  V  A  2 

- v  -  =  -  >  J,  Jb  ~  — -v  >  J,  therefore,  the  a's  can  be 

B2  J2  100  B2 


neglected.  Let 


2  2 
A  2  -2  A  -2  mr?  A 

b  =  v  +4Jc  JM,  and  c  =  2  —  v  [  JN-JD~k  JH~  {2JT~JA)  — =— ]  =2  — v-f. 


B 


A 


B 


such  that  now  we  can  rewrite  inequality  A. 54  as 
2  _  _  (f24v)2  _  4  B2k~2J„  j=2 


A".. 2  -2  T  .  <f2|V)2.2  .  ,  f. 

2 V  AT  4  a  _ ,  2.. 2  2  (A. 55) 


B 


aA  v 


1  4B2k~2JM 

Recalling  the  stability  factor,  S,  —  = -  allows  the 


following  steps: 


A  v 


/  -r  t  ,  -2  _  T  md  ,  _  _ 
a-  (  JD+JN+k  JH  JA~ ~S1 


-2 


md 


a+f-  (  J^+J^j+k  J  u  J  a — - — +  ( Jj\j  J  p,  k  J  (2Jrp  J  %) 


md‘ 


H  A 


N  UD 


*  fj  \^u  ijr>  u 
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,2  -2  md 2 

a-f={~Jci+JN+k  JH~  ~ mS3 


These  terms  S-j_,  S2 ,  S3  are  very  important  to  the  rest  of  the 
discussion  of  stability.  The  final  criteria  for  stability  is: 


1  <  Elh. 

s  si 


S1  >  0 


(A. 56; 


.  -2  md‘ 

where  Sx  =  JN+k  JH~JD~J a~~X 


52  =  2JN  2JD_2JT  ~a~> 

53  =  2k~2JH+{2JT~2JA)^~ 


(A. 57 ) 


Using  the  relationships  between  the  S  terms  and  f,  the 
stability  criteria  can  be  rewritten. 


S2S3  (S-L  +  f)  (Sx-f)  s2_f2 


s1>0 


(A. 58) 


-|J  [  1y fJ^S)  ±y/]o+Tc]  dp 


Then,  A«e 


-2  ,  A  2  „  ,  -2 


where  b=4Jc  Jw- — — v  =4ic  JM(1- 


a2v2 


,M  —v  —  _2 


)  =  4ic 


4S  Jc 


C=2^Vf  =  2^4ic  2Jm  JS-f 


and  because  S>0,  it  follows  that  S>1  from  A.  58.  This  is  the 

"classical  stability  condition."  This  is  not  sufficient  though. 

The  general  motion  for  . r  , 

-i  I  v  dp 
l  =  Ule  J  , 


with  the  two  solutions  for  A  represents  the  superposition  of  the 
two  cyclic  motions,  each  decreasing  in  amplitude.  This  is  the 
epicyclic  motion  caused  by  the  coning  and  spinning  simultaneously. 
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For  S  values  around  1.5,  the  epicyclic  yawing  motion  will  have  many 
maxima.  The  successive  values  of  the  yaw  at  maxima  may  increase 
without  bound. 

If  the  classical  condition  for  stability  is  not  met  (S<1)  , 
then  the  final  step  of  the  inequality  A. 58  can  not  be  satisfied. 
The  yaw  will  grow  steadily  and  the  motion  is  similar  to  the  falling 
motion  of  a  top  as  the  spin  becomes  too  small. 

It  would  be  difficult  to  meet  the  condition  S^CO  since  the 
terms  JN,  10 JD.  If  JH  <  0,  then  this  condition  could  be  met, 
but  this  is  unlikely  to  happen  since  it  would  require  a  shift  of 
the  center  of  mass,  which  for  our  purposes  is  not  practical. 

Finally,  the  stability  found  in  the  product  of  the  S 2,  S3 
terms  may  be  sufficient  to  cause  instability.  Recall  that: 


1  F  F  °2  f 

—  <  l- —  for  S-,  >0,  — -  could  become  >  l.  — —  =1  +  -— , 
S  2  1  s2  Sl  Sl 


This  would  happen  if  either  of  the  S2,  S3  terms  become  negative. 
But,  because  S2+S2~2S^,  S2+S3>0,  always.  Depending  on  the  sign  of 
the  Magnus  moment  coefficient,  Jj, ,  a  change  in  magnitude  could 
cause  either  S  term  to  become  negative  and  result  in  an  unstable 
projectile. 
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The  required  heat  is  best  deposited  in  a  time,  tf  <  tD.  The 
average  temperature  to  which  the  material  must  be  raised  is 
estimated  to  be: 


-  '^M+r^s 

T  =  =  2700  C. 


(B .  2 ) 


This  gives  a  required  f luence  of : 

F0tf  =  p£[LM+  (T-Tq)  Cs  =  4700  J/ cm2  .  (B.3) 

For  continuous  laser  illumination  with  a  sufficiently  large 
burn  spot  (radius  30  cm)  on  a  projectile  rotating  with  a  period  of 
4  milliseconds,  an  equal  distribution  of  the  f luence  is  achieved. 

The  absorptance  of  steel  is  approximately  30%  (Schriempf, 
1974)  .  This  means  that  for  a  laser  intensity,  IL,  the  f luence 
becomes : 

F0t  f  =  tfO  .3-1 L  =  4.1  KJ/ cm2  .  (b  .  4 ) 


This  gives  a  required  laser  energy  of 


nr‘ 


„  4700  ,  2 

=  2* - .3  J/cm  , 


(B .  5) 


where  the  factor  2  takes  into  account  the  energy  required  on  both 
sides  of  the  projectile  and  the  factor  0.3  accounts  for  the 
reflectance.  The  result  of  solving  equation  B.5  for  the  energy  is 
nearly  90  mega  joules  required.  For  a  tf  <  0.45  sec,  this 
corresponds  to  a  laser  power  of  200  megawatts. 

If  the  burn  spot  is  reduced  to  an  8  centimeter  radius  and  it 
is  assumed  the  laser  energy  is  moved  over  the  projectile  body,  then 
substituting  values  into  equation  B.5  results  in  an  energy 
requirement  of  only  6  mega joules  and  a  power  of  14  megawatts. 
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